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2 DIFFERENTIAL EQUATIONS 


Differential Equations 


A differential equation is an equation involving an unknown function and 
its derivatives. 


Example 1.1: The following are differential equations involving the un- 
known function y. 


Dea (1.1) 
dx 
2 2 
oh 4( 2) a1 (1.2) 
x Ix 
d° ee a 
473 + (sinx) 5 + 5ay=0 (1.3) 
3 
d’y (#) 12) 
—> | +3y, =] ty} —] =5 1.4 
(3) age 4 dx ae) 
dy dy 
pacer ite, Ian ae 9) 1.5 
or? ax? a>) 


A differential equation is an ordinary differential equation if the unknown 
function depends on only one independent variable. If the unknown func- 
tion depends on two or more independent variables, the differential equa- 
tion is a partial differential equation. In this book we will be concerned 
solely with ordinary differential equations. 


Example 1.2: Equations 1.1 through 1.4 are examples of ordinary differ- 
ential equations, since the unknown function y depends solely on the vari- 
able x. Equation 1.5 is a partial differential equation, since y depends on 
both the independent variables t and x. 
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he Note! 


The order of a differential equation is the order of 
the highest derivative appearing in the equation. 


Example 1.3: Equation 1.1 is a first-order differential equation; 1.2, 1.4, 
and 1.5 are second-order differential equations. (Note in 1.4 that the or- 
der of the highest derivative appearing in the equation is two.) Equation 
1.3 is a third-order differential equation. 


Notation 


The expressions y’,y”,y”,y™,...,y are often used to represent, re- 
spectively, the first, second, third, fourth, . . ., nth derivatives of y with re- 
spect to the independent variable under consideration. Thus, y” repre- 


sents d? y/ dx” if the independent variable is x, but represents d? yl dp” 


if the independent variable is p. Observe that parenthesis are used in y”’ 
to distinguish it from the nth power, y’. If the independent variable is 
time, usually denoted by ¢, primes are often replaced by dots. Thus, 


y, ¥, and ¥ represent, dy/ dt, d 2 y/ dt’, and d 3y i dt’, respectively. 


Solutions 


A solution of a differential equation in the unknown function y and the 
independent variable x on the interval # is a function y(x) that satisfies 
the differential equation identically for all x in §. 


Example 1.4: Is y(x)=c, sin2x +c, cos2x, where c, and c, are arbi- 
trary constants, a solution of y”+4y=0? 


Differentiating y, we find y’ = 2c,cos2x — 2c, sin2x and y" = 
—4c, sin2x — 4c, cos2x . Hence, 
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y’+4y = (—4c,sin2x—4c, cos2x)+4(c, sin2x +c, cos2x) 
= (—4c, +4c,)sin2x + (—4c, + 4c, )cos2x 
= 0 


Thus, y=c,sin2x+c, cos2x satisfies the differential equation for all 
values of x and is a solution on the interval (—°0,00). 


Example 1.5: Determine whether y= x* —1 is a solution of (y’)* +? 
=-l. 


Note that the left side of the differential equation must be nonnegative for 
every real function y(x) and any x, since it is the sum of terms raised to 
the second and fourth powers, while the right side of the equation is neg- 
ative. Since no function y(x) will satisfy this equation, the given differ- 
ential equation has no solutions. 

We see that some differential equations have infinitely many solu- 
tions (Example 1.4), whereas other differential equations have no solu- 
tions (Example 1.5). It is also possible that a differential equation has ex- 
actly one solution. Consider (y’)4 + y = 0, which for reasons identical 
to those given in Example 1.5 has only one solution y=0. 


You Need to Know J 


A particular solution of a differential equation is any 
one solution. The general solution of a differential 
equation is the set of all solutions. 


Example 1.6: The general solution to the differential equation in Ex- 
ample 1.4 can be shown to be (see Chapters Four and Five) y= 
c, sin2x +c, cos2x. That is, every particular solution of the differential 
equation has this general form. A few particular solutions are: (a) y= 
5sin2x—3cos2x (choose c;=5 and c, =—3), (b) y=sin2x (choose 
c, =1 and c, =0), and (c) y=0 (choose c, =c, =0). 
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The general solution of a differential equation cannot always be ex- 
pressed by a single formula. As an example consider the differential equa- 


tion y’+ y = 0, which has two particular solutions y=1/x and y=0. 


Initial-Value and Boundary-Value Problems 


A differential equation along with subsidiary 
conditions on the unknown function and its de- 
rivatives, all given at the same value of the in- 
dependent variable, constitutes an initial-value 
problem. The subsidiary conditions are initial 
conditions. If the subsidiary conditions are giv- 
en at more than one value of the independent 
variable, the problem is a boundary-value prob- 
lem and the conditions are boundary conditions. 


Example 1.7: The problem y” +2y’ =e*;y(z)=1,y’(a) =2 is an initial 
value problem, because the two subsidiary conditions are both given at 
x=a.The problem y”+2y’ = e*;y(0)=1, y(1) =1 is a boundary-value 
problem, because the two subsidiary conditions are given at x = 0 and 
x=1. 

A solution to an initial-value or boundary-value problem is a func- 
tion y(x) that both solves the differential equation and satisfies all given 
subsidiary conditions. 


Standard and Differential Forms 


Standard form for a first-order differential equation in the unknown func- 
tion y(x) is 


y= f(xy) (1.6) 


where the derivative y’ appears only on the left side of 1.6. Many, but 
not all, first-order differential equations can be written in standard form 
by algebraically solving for y’ and then setting f(x,y) equal to the right 
side of the resulting equation. 
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The right side of 1.6 can always be written as a quotient of two oth- 
er functions M(x,y) and —N(x,y). Then 1.6 becomes dy/dx = M(x, y)/ 
—N(x,y), which is equivalent to the differential form 


M(x, y)dx + N(x, y)dy =0 (1.7) 


Linear Equations 


Consider a differential equation in standard form 1.6. If f(x,y) can be writ- 
tenas f(x, y)=—p(x)y+q(x) (that is, as a function of x times y, plus an- 
other function of x), the differential equation is linear. First-order linear 
differential equations can always be expressed as 


y+ p(x)y = q(x) (1.8) 
Linear equations are solved in Chapter Two. 
Bernoulli Equations 
A Bernoulli differential equation is an equation of the form 
y+ pQx)y=qa)y" (1.9) 


where n denotes a real number. When n = | or n=0, a Bernoulli equation 
reduces to a linear equation. Bernoulli equations are solved in Chapter 
Two. 


Homogeneous Equations 
A differential equation in standard form (1.6) is homogeneous if 
f(tx,ty) = f(x,y) (1.10) 


for every real number ¢. Homogeneous equations are solved in Chapter 
Two. 
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he Note! 


In the general framework of differential equations, 
the word “homogeneous” has an entirely different 
meaning (See Chapter Four). Only in the context 
of first-order differential equations does “homoge- 
neous” have the meaning defined above. 


Separable Equations 
Consider a differential equation in differential form (1.7). If M(x,y) =A(@) 
(a function only of x) and N(x,y) = B(y) (a function only of y), the differ- 


ential equation is separable, or has its variables separated. Separable 
equations are solved in Chapter Two. 


Exact Equations 
A differential equation in differential form (1.7) is exact if 


OM(x,y) _ AN(x,y) 
oy ox 


(1.11) 


Exact equations are solved in Chapter Two (where a more precise defini- 
tion of exactness is given). 


Chapter 2 
SOLUTIONS OF 
FIRST-ORDER 


DIFFERENTIAL 
EQUATIONS 


IN THIS CHAPTER: 


VY Separable Equations 
Homogeneous Equations 
Exact Equations 

Linear Equations 
Bernoulli Equations 

Y Solved Problems 


V 
Vv 
v 
Vv 


Separable Equations 
General Solution 


The solution to the first-order separable differential equation (see Chap- 
ter One). 


A(x)dx + B(y)dy =0 (2.1) 
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J A@ax de J B(y)dy =c (2.2) 


where c represents an arbitrary constant. 

(See Problem 2.1) 

The integrals obtained in Equation 2.2 may be, for all practical pur- 
poses, impossible to evaluate. In such case, numerical techniques (see 
Chapter 14) are used to obtain an approximate solution. Even if the indi- 
cated integrations in 2.2 can be performed, it may not be algebraically 
possible to solve for y explicitly in terms of x. In that case, the solution is 
left in implicit form. 


Solutions to the Initial-Value Problem 
The solution to the initial-value problem 

A(x)dx + B(y)dy=0;  y(x9) = yo (2.3) 
can be obtained, as usual, by first using Equation 2.2 to solve the differ- 
ential equation and then applying the initial condition directly to evalu- 


ate c. 
Alternatively, the solution to Equation 2.3 can be obtained from 


x y 
J A(s)ds +f B(t)dt = 0 (2.4) 
Xo Yo 


where s and ¢ are variables of integration. 


Homogeneous Equations 

The homogeneous differential equation 
d 
= f(xy) (2.5) 
dx 


having the property f(tx, ty) = f(x, y) (see Chapter One) can be trans- 
formed into a separable equation by making the substitution 
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y=xv (2.6) 
along with its corresponding derivative 
= =vtx— (2.7) 
The resulting equation in the variables v and x is solved as a separable 
differential equation; the required solution to Equation 2.5 is obtained by 
back substitution. 


Alternatively, the solution to 2.5 can be obtained by rewriting the dif- 
ferential equation as 


ox 2.8 
dy f(x,y) = 
and then substituting 
x= yu (2.9) 
and the corresponding derivative 
at =ut ye (2.10) 
dy dy 


into Equation 2.8. After simplifying, the resulting differential equation 
will be one with variables (this time, u and y) separable. 

Ordinarily, it is immaterial which method of solution is used. Occa- 
sionally, however, one of the substitutions 2.6 or 2.9 is definitely superi- 
or to the other one. In such cases, the better substitution is usually appar- 
ent from the form of the differential equation itself. 

(See Problem 2.2) 


Exact Equations 
Defining Properties 
A differential equation 


M(x, y)dx + N(x, y)dy = 0 (2.11) 
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is exact if there exists a function g(x, y) such that 


dg(x, y) = M(x, y)dx + N(x, y)dy (2.12) 


ke Note! 


Test for exactness: lf M(x,y) and N(x,y) are con- 
tinuous functions and have continuous first partial 
derivatives on some rectangle of the xy-plane, then 
Equation 2.11 is exact if and only if 
OM(x,y) _ ON(x, y) 
oy ox 


(2.13) 


Method of Solution 


To solve Equation 2.11, assuming that it is exact, first solve the equations 


OY) _ Mx, y) (2.14) 
Ox 

Og(x,y) _ 

ae = N(x, y) (2.15) 


for g(x, y). The solution to 2.11 is then given implicitly by 
g(x, y)=c (2.16) 


where c represents an arbitrary constant. 

Equation 2.16 is immediate from Equations 2.11 and 2.12. If 2.12 is 
substituted into 2.11, we obtain dg(x, y(x)) = 0. Integrating this equation 
(note that we can write 0 as 0 dx), we have Jdg(x,.yOo) = JOax, which, 
in turn, implies 2.16. 


Integrating Factors 
In general, Equation 2.11 is not exact. Occasionally, it is possible to trans- 


form 2.11 into an exact differential equation by a judicious multiplica- 
tion. A function J(x, y) is an integrating factor for 2.11 if the equation 
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I(x, y)[M(x, y)dx + N(x, y)dy] = 0 (2.17) 
is exact. A solution to 2.11 is obtained by solving the exact differential 


equation defined by 2.17. Some of the more common integrating factors 
are displayed in Table 2.1 and the conditions that follow: 


(#4 - wy) = g(x), a function of x alone, then 
y x 
(x,y) =e)" (2.18) 
If mike - wy) = h(y), a function of y alone, then 
y x 
(=e O* (2.19) 
If M = yf(xy) and N = xg(xy), then 
I(x,y)=——— (2.20) 


xM —yN 


In general, integrating factors are difficult to uncover. If a differential 
equation does not have one of the forms given above, then a search for 
an integrating factor likely will not be successful, and other methods of 
solution are recommended. 

(See Problems 2.3—2.6) 


Linear Equations 


Method of Solution 


A first-order linear differential equation has the form (see Chapter One) 


y’ + p(x)y = q(x) (2.21) 


An integrating factor for Equation 2.21 is 


eo 


I(x)=e (2.22) 
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Group of terms 


Table 


Integrating factor /(x, y) 


2.1 


Exact differential dy(x, v) 


yvdx —xdy 


+- 


xdy—yd 
xdy ~ydx _ fy 


2 


ydx —xdy 


ydx —xdy 


vdx —xdy 


ydx +xdy 


d(In xy) 


ydx +xdy 


ydx +x dy _ 
(xy)" 


lier 


vdytxdx 


; =al5 in (e+y)| 


x+y 2 


ydy+xdx 


ydy +xdx 


a: me im ae 
Ge + yy" 


= d| - Sens 
ce “Deetyy! 


ay dx + bx dy 
(a, b constants) 


x“ ly? "ay dx + bx dy) = d(x“y’) 


which depends only on x and is independent of y. When both sides of 2.21 
are multiplied by /(x), the resulting equation 


I(x)y’ + p(x)I(x 


)y = I(x)q(x) (2.23) 


is exact. This equation can be solved by the method described previous- 
ly. A simpler procedure is to rewrite 2.23 as 


d(yD) _ 
dx 


Iq(x) 
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integrate both sides of this last equation with respect to x, and then solve 
the resulting equation for y. The general solution for Equation 2.21 is 


_ J I(x)q(x)dx +e 
I(x) 


where c is the constant of integration. 
(See Problem 2.7) 


Bernoulli Equations 
A Bernoulli differential equation has the form 
y+ p@x)y=qia)y" (2.24) 
where 7 is a real number. The substitution 
Zaye (2.25) 
transforms 2.24 into a linear differential equation in the unknown func- 


tion z(x). 
(See Problem 2.8) 


Solved Problems 
2 
Solved Problem 2.1 Solve & = * +? 
dx y 
This equation may be rewritten in the differential form 
(x? +2)dx — ydy =0 
which is separable with A(x) = x? + 2 and B(y) =—y. Its solution is 
Joe +2)dx- | ydy=c 
or 


1 342,-1y =c 
3 2 
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Solving for y, we obtain the solution in implicit form as 
2 
y= a xe+4x+k 


with k =—2c. Solving for y explicitly, we obtain the two solutions 


y= Seder and y=- Seder 


+x 


Solved Problem 2.2 Solve y’ =~ 


This differential equation is not separable. Instead it has the form 
y’ = f(x,y), with 


yx 


f(x y= 


where 


ty+tx  ty+x)_ y+x 
tx tx x 


f(tx, ty) = = f(y) 


so it is homogeneous. Substituting Equations 2.6 and 2.7 into the equa- 
tion, we obtain 


dv xv+x 
vtx— = —— 
dx x 


which can be algebraically simplified to 


This last equation is separable; its solution is 
1 
—dx-|dv=c 
J—ae-] 
which, when evaluated, yields v =In|x|—c, or 
v=InIkx| (2.26) 
where we have set c = —In|k |and have noted that In| x«1+In|k|=In| xk |. 


Finally, substituting v = y/x back into 2.26, we obtain the solution to the 
given differential equation as y= xIn|kx|. 
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Solved Problem 2.3 Solve 2xydx+(1+x7)dy=0. 


This equation has the form of Equation 2.11 with M(x, y) = 2xy and 
N(x, y) = 1+ x2. Since dM / dy =0N/ 0x = 2x, the differential equation 
is exact. Because this equation is exact, we now determine a function 
g(x, y) that satisfies Equations 2.14 and 2.15. Substituting M(x, y) = 2xy 
into 2.14, we obtain dg / dx = 2xy. Integrating both sides of this equation 
with respect to x, we find 


dg 
J oe dx =| 2xydx 
or 
g(x,y) =x y+ hy) (2.27) 
Note that when integrating with respect to x, the constant (with respect to 
x) of integration can depend on y. 

We now determine hA(y). Differentiating 2.27 with respect to y, we 
obtain dg / dy = x” +h’(y) Substituting this equation along with N(x, y) = 
1 +x? into 2.15, we have 

2 , = 2 , a 
x +h'(y)=1+x° or h’(y)=1 


Integrating this last equation with respect to y, we obtain h(y) = y+ c, 
(c, = constant). Substituting this expression into 2.27 yields 


gi, y)=x*y+y+e, 


The solution to the differential equation, which is given implicitly by 
2.16 as g(x, y) =c, is 


xyty=c, (c,=c-c¢,) 
Solving for y explicitly, we obtain the solution as y = c,/(x? + 1). 


Solved Problem 2.4 Determine whether the differential equation ydx — 
xdy = 0 is exact. 
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This equation has the form of Equation 2.11 with M (x, y)=y and N (x, y) 
=—x. Here 


as and ON, 54 


oy ox 
which are not equal, so the differential equation is not exact. 


Solved Problem 2.5 Determine whether —1/x? is an integrating factor 
for the differential equation ydx — xdy = 0. 


It was shown in Problem 2.4 that the differential equation is not exact. 
Multiplying it by —1/x?, we obtain 


S (ydx —xdy) =0 or a+ dy=0 (2.28) 
x x xX 


Equation 2.28 has the form of Equation 2.11 with M (x, y) = —y/x? and 
N (x, y) = 1/x. Now 


oM_O (+)-5 _ oO (=)-= 

dy dykx?) x? axlx) ax 
so 2.28 is exact, which implies that —1/x? is an integrating factor for the 
original differential equation. 


Solved Problem 2.6 Solve ydx — xdy =0. 


Using the results of Problem 2.5, we can rewrite the given differential 
equation as 
dy — yd 
xdy - ag) 


x 


which is exact. Equation 2.28 can be solved using the steps described in 
Equations 2.14 through 2.16. 

Alternatively, we note from Table 2.1 that 2.28 can be rewritten as 
d (y/x) = 0. Hence, by direct integration, we have y/x =c, or y = xc, as 
the solution. 


Solved Problem 2.7 Solve y’+(4/x)y= x". 
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The differential equation has the form of Equation 2.21, with p(x) = 4/x 
and q(x) = x‘, and is linear. Here 


4 
J p@oax = Jrax =4InIx!l=Inx* 
x 
so 2.22 becomes 
ese OP ae ee (2.29) 


Multiplying the differential equation by the integrating factor de- 
fined by 2.29, we obtain 


; d 
x+y +4x3y =x® or —(yx*) =x 
dx 
Integrating both sides of this last equation with respect to x, we obtain 


4 


x Sees or sig is 
coer ae a) 


Solved Problem 2.8 Solve y’+xy = xy’. 


This equation is not linear. It is, however, a Bernoulli differential equa- 
tion having the form of Equation 2.24 with p(x) = g(x) =x, andn=2. We 
make the substitution suggested by 2.25, namely z = y!-? = y"!, from 
which follow 
y= 2 and y'=- aoe 
z z 


Substituting these equations into the differential equation, we obtain 


This last equation is linear for the unknown function z(x). It has the form 
of Equation 2.21 with y replaced by z and p(x) = g(x) = —x. The integrat- 
ing factor is 


I(x)= gone =e? 


Multiplying the differential equation by /(x), we obtain 
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-2/2 az a) _x? 
e w12 ©* _ ye x 27 _ xe x°/2 
dx 


or 
d re _2 
“(ze x 2) = xe x72 
dx 

Upon integrating both sides of this last equation, we have 


2 2 
ze ae Pa4e 


whereupon 


12 


z(x)= ce* 2 44 


The solution of the original differential equation is then 
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VSN AN 


Growth and Decay Problems 


Let M(t) denote the amount of substance (or popula- 


tion) that is either growing or decaying. If we assume 
that dN/dt, the time rate of change of this amount of 
substance, is proportional to the amount of substance 


present, then dN/dt = kN, or 
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ans kN=0 (3.1) 
dt 
where k is the constant of proportionality. 

We are assuming that M(f) is a differentiable, hence continuous, func- 
tion of time. For population problems, where N(‘) is actually discrete and 
integer-valued, this assumption is incorrect. Nonetheless, 3.1 still pro- 
vides a good approximation to the physical laws governing such a sys- 
tem. 


Temperature Problems 


Newton’s law of cooling, which is equally applicable to heating, states 
that the time rate of change of the temperature of a body is proportional 
to the temperature difference between the body and its surrounding medi- 
um. Let T denote the temperature of the body and let 7, denote the tem- 
perature of the surrounding medium. Then the time rate of change of the 
temperature of the body is dT/dt, and Newton’s law of cooling can be for- 
mulated as dT/dt = —k(T — T_,), or as 


a KT, (3.2) 
dt 


where k is a positive constant of proportionality. Once k is chosen posi- 
tive, the minus sign is required in Newton’s law to make dT/dt negative 
in a cooling process, when 7 is greater than T_, and positive in a heating 
process, when 7 is less than T_. 


Falling Body Problems 


Consider a vertically falling body of mass m that is being influenced only 
by gravity g and an air resistance that is proportional to the velocity of 
the body. Assume that both gravity and mass remain constant and, for 
convenience, choose the downward direction as the positive direction. 
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You Need to Know J 


Newton’s second law of motion: The net force 
acting on a body is equal to the time rate of change 
of the momentum of the body; or, for constant 
mass, 


F=m— (3.3) 


where F is the net force on the body and v is the 
velocity of the body, both at time t. 


For the problem at hand, there are two forces acting on the body: the force 
due to gravity given by the weight w of the body, which equals mg, and 
the force due to air resistance given by —kv, where k = 0 is a constant of 
proportionality. The minus sign is required because this force opposes the 
velocity; that is, it acts in the upward, or negative, direction (see Figure 
3-1). The net force F on the body is, therefore, F = mg — kv. Substituting 
this result into 3.3, we obtain 


dv 
_kbw=m— 
oo mat 
or 
OTB ey (3.4) 
dt m 


as the equation of motion for the body. 
If air resistance is negligible or nonexistent, then k = 0 and 3.4 sim- 
plifies to 
av 
dt 
When k > 0, the limiting velocity v, is defined by 


& (3.5) 


v= oe (3.6) 
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kv 
Falling body 
v mg 


Ground 
Y 
Y Uj 
Positive x-direction 


Figure 3-1 


Caution: Equations 3.4, 3.5, and 3.6 are valid only if the given con- 
ditions are satisfied. These equations are not valid if, for example, air re- 
sistance is not proportional to velocity but to the velocity squared, or if 
the upward direction is taken to be the positive direction. 


Dilution Problems 


Consider a tank which initially holds V, gal of brine that contains a Ib of 
salt. Another solution, containing b lb of salt per gallon, is poured into the 
tank at the rate of e gal/min while simultaneously, the well-stirred solu- 
tion leaves the tank at the rate of f gal/min (Figure 3-2). The problem is 
to find the amount of salt in the tank at any time t. 

Let Q denote the amount (in pounds) of salt in the 
tank at any time. The time rate of change of Q, dQ/dt, 
equals the rate at which salt enters the tank minus the rate 
at which salt leaves the tank. Salt enters the tank at the 
rate of be lb/min. To determine the rate at which salt 
leaves the tank, we first calculate the volume of brine 
in the tank at any time t, which is the initial volume V, 
plus the volume of brine added et minus the volume of brine removed ft. 
Thus, the volume of brine at any time is 
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VY) t+et—ft (3.7) 


The concentration of salt in the tank at any time is Q/ (Vj + et — ft), from 
which it follows that salt leaves the tank at the rate of 


(25) Ib/min 
Vo +et— ft 


40 _ 5p (aes) 
dt VY tet—ft 


i ns 
dt Vy+(e-f)t 


Thus, 


or 


be (3.8) 


Figure 3-2 
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R 


Figure 3-3 


Electrical Circuits 


The basic equation governing the amount of current J (in amperes) in a 
simple RL circuit (see Figure 3-3) consisting of a resistance R (in ohms), 
an inductor L (in henries), and an electromotive force (abbreviated emf) 
E (in volts) is 

dl fe R i= £ 3.9 

dt LL ee 
For an RC circuit consisting of a resistance, a capacitance C (in farads), 
an emf, and no inductance (Figure 3-4), the equation governing the 
amount of electrical charge g (in coulombs) on the capacitor is 

i Ease (3.10) 


dt RC! R 


Figure 3-4 


26 DIFFERENTIAL EQUATIONS 
The relationship between g and J is 
(3.11) 


For more complex circuits see Chapter Seven. 


Orthogonal Trajectories 
Consider a one-parameter family of curves in the xy-plane defined by 
F(x, y, c) =0 (3.12) 


where c denotes the parameter. The problem is to find another one-pa- 
rameter family of curves, called the orthogonal trajectories of the fami- 
ly of curves in 3.12 and given analytically by 


G(x, y, k) =0 (3.13) 


such that every curve in this new family 3.13 intersects at right angles 
every curve in the original family 3.12. 

We first implicitly differentiate 3.12 with respect to x, then eliminate 
c between this derived equation and 3.12. This gives an equation con- 
necting x, y, and y’, which we solve for y’ to obtain a differential equa- 
tion of the form 


d 
& = f(x,y) (3.14) 
dx 
The orthogonal trajectories of 3.12 are the solutions of 
d___it (3.15) 
dx f(x,y) 


For many families of curves, one cannot explicitly solve for dy/dx 
and obtain a differential equation of the form 3.14. We do not consider 
such curves in this book. 


Solved Problems 


Solved Problem 3.1 A bacteria culture is known to grow at a rate pro- 
portional to the amount present. After one hour, 1000 strands of the bac- 
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teria are observed in the culture; and after four hours, 3000 strands. Find 
(a) an expression for the approximate number of strands of the bacteria 
present in the culture at any time ¢ and (b) the approximate number of 
strands of the bacteria originally in the culture. 


(a) Let N(t) denote the number of bacteria strands in the culture at 
time ¢. From Equation 3.1, dN/dt — kN = 0, which is both linear and sep- 
arable. Its solution is 


N(t) = ce (3.16) 
At t= 1, N= 1000; hence, 

1000 = ce* (3.17) 
At t= 4, N = 3000; hence. 

3000 = ce** (3.18) 


Solving 3.17 and 3.18 for k and c, we find 
k= zin3 ~ 0.3662 and c= 1000e~* = 693 
Substituting these values of k and c into 3.16, we obtain 
N(t) = 693e93667" (3.19) 


(b) We require N at t= 0. Substituting t= 0 into 3.19, we obtain N(0) 
= 693 e0-3662)) — 693, 


Solved Problem 3.2 A tank initially holds 100 gal of a brine solution 
containing 20 Ib of salt. At t= 0, fresh water is poured into the tank at the 
rate of 5 gal/min, while the well-stirred mixture leaves the tank at the 
same rate. Find the amount of salt in the tank at any time t. 


Here, V, = 100, a = 20, b = 0, and e = f= 5. Equation 3.8 becomes 


dQ 1 
—~=+—Q=0 
dt 30 2 
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The solution of this linear equation is 

O = ce/20 (3.20) 
At t=0, we are given that Q = a= 20. Substituting these values into 3.20, 


we find that c = 20, so that 3.20 can be rewritten as O = 20e~*/?°. Note 
that as t > o, Q > 0 as it should, since only fresh water is being added. 
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An nth-order linear differential equation has the form 


b, (xy + by yO? +2 +B Dy’ + boy = s(x) (4.1) 


where g(x) and the coefficients b(x) (j=0,1,2,...,) de- 
pend solely on the variable x. In other words, they do 
not depend on y or any derivative of y. 
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If g(x) = 0, then Equation 4.1 is homogeneous; if not, 4.1 is nonho- 
mogeneous. A linear differential equation has constant coefficients if all 
the coefficients b(x) in 4.1 are constants; if one or more of these coeffi- 
cients is not constant, 4.1 has variable coefficients. 


Theorem 4.1. Consider the initial-value problem given by the linear dif- 
ferential equation 4.1 and the n initial conditions 


WX=Co, WA)=Q, 


; ae (4.2) 
y (x9) = Copeeeny D(x) =Cy-y 


If g(x) and b(x) (j = 0,1,2,..., 2) are continuous in some interval # con- 
taining x, and if b (x) #0 in #, then the initial-value problem given by 
4.1 and 4.2 has a unique (only one) solution defined throughout #. 


When the conditions on b(x) in Theorem 4.1 hold, we can divide 
Equation 4.1 by b (x) to get 


y” +4, (x) yr tes +d5(x)y” +a, (x)y’ + dg (x)y = (x) (4.3) 


where a(x) = b (xb, (x) (j = 0,1,2,..., 2 — 1) and O(x) = g(x)/b, (x). 
Let us define the differential operator L(y) by 


L(y) = yy +a, (x)yO? +++ ay (x)y” +a, (xy tay(x)y (4.4) 


where a,(x) (i= 0,1,2,..., 2 — 1) is continuous on some interval of interest. 
Then 4.3 can be rewritten as 


L(y) = @) (4.5) 


and, in particular, a linear homogeneous differential equation can be ex- 
pressed as 


LQ) =0 (4.6) 
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Linearly Independent Solutions 


A set of functions {y, (x), y>(%),-..,y,,(%)} 1s linearly dependent ona Sx Sb 
if there exist constants c,,c,,...,c,, not all zero, such that 
CY (4) + Cp 9 (x) +++ +¢,y, (4) = 0 (4.7) 


onasx<b. 


Example 4.1: The set {x,5x,1,sin x}is linearly dependent on [—1,1] since 
there exist constants c, = —5, c, = 1, c, = 0, and c, = 0, not all zero, such 
that 4.7 is satisfied. In particular, 


—5-x+1-5x+0-14+0-sinx=0 


Note that c, = c,=--- c, =O is a set of constants that always satis- 
fies 4.7. A set of functions is linearly dependent if there exists another set 
of constants, not all zero, that also satisfies 4.7. If the only solution to 4.7 
iS c, = c,=-+-c, = 0, then the set of functions {y,(x), y,(4),....¥,,(%)} is 
linearly independent onas x <b. 


Theorem 4.2. The nth-order linear homogeneous differential equation 


L(y) = 0 always has n linearly independent solutions. If y, (x),y,(4),.-..y,,%) 
represent these solutions, then the general solution of L(y) = 0 is 


y(X) = CY, (0) + Cry) +++ Cy, Or) (4.8) 
where c,,c 


pro 


The Wronskian 


c,, denote arbitrary constants. 


The Wronskian of a set of functions {z,(x), z,(),..., Z,(«)} on the interval 
a<x<b, having the property that each function possesses n — J deriva- 
tives on this interval, is the determinant 


ra Zt ky 
, , i 
zi Bow 2g 
Ww aw aw 
W(Z.Zo.0Z I=] y aN Cn 


zn) zr) nae ga) 
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Theorem 4.3. If the Wronskian of a set of n functions defined on the in- 
terval a <x <b is nonzero for at least one point in this interval, then the 
set of functions is linearly independent there. If the Wronskian is identi- 
cally zero on this interval and if each of the functions is a solution to the 
same linear differential equation, then the set of functions is linearly de- 
pendent. 


Caution: Theorem 4.3 is silent when the Wronskian is identically 
zero and the functions are not known to be solutions of the same linear 
differential equation. In this case, one must test directly whether Equa- 
tion 4.7 is satisfied. 


Nonhomogeneous Equations 

Let Ms denote any particular solution of Equation 4.5 (see Chapter One) 
and let y, (henceforth called the homogeneous or complementary solu- 
tion) represent the general solution of the associated homogeneous equa- 


tion L(y) = 0. 


Theorem 4.4. The general solution to L(y) = ¢(x) is 


y=y,tY, (4.9) 


Don’t Forget 


The general solution is the sum of 
the homogeneous and particular so- 
lutions. 
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The Characteristic Equation 


Second-Order Equations 
Corresponding to the differential equation 

y’+ay +ajy=0 (5.1) 
in which a, and ay are constants, is the algebraic equation 

MV +aA+ay =0 (5.2) 
which is obtained by substituting y = e** (assuming x to be the indepen- 
dent variable) into Equation 5.1 and simplifying. Note that Equation 5.2 
can be easily obtained by replacing y”, y’, and y by A”, 4}, and, 2° = 1, 


respectively. Equation 5.2 is called the characteristic equation of 5.1. 


Example 5.1 The characteristic equation of y” +3y’—4y =0 isA*+3A 
—4=0; the characteristic equation of y”—2y’+y=0 isA?—22+1=0. 


Characteristic equations for differential equations having dependent 
variables other than y are obtained analogously, by replacing the jth de- 
rivative of the dependent variable by A/ (j = 0,1,2). 
nth-Order Equations 
Similarly, the characteristic equation of the differential equation 

Og ja) Aas oc wat - 
y Any +-+ay +dajy=0 (5.3) 
with constant coefficients a; (j=0,1,...,n — 1) is 
A" +4, An +++ ajA+dy =0 (5.4) 
The characteristic equation 5.4 is obtained from 5.3 by replacing y by 
Aj (j = 0,1,..., ). Characteristic equations for differential equations hav- 
ing dependent variables other than y are obtained analogously, by 


replacing the jth derivative of the dependent variable by 4/ (j = 0,1...., 
n). 
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mt 


Example 5.2. The characteristic equation of y —3y’” +2y”—y=0 is 


A4 — 343 + 2A? — 1 =0. The characteristic equation of 


A =347°45A=7=0 


kK Caution! 


Characteristic equations are only defined for linear 
homogeneous differential equations with constant 
coefficients. 


General Solution for Second-Order Equations 
The characteristic equation (5.2) can be factored into 


(A-A )A-Aq) = 0 (5.5) 


wx Note! 


The roots of the characteristic polynomial deter- 
mine the solution of the differential equation. 


There are three cases to consider. 


Case 1. 4, and 4, both real and distinct. Two linearly indepen- 
dent solutions are e”'* and e”*, and the general solution is (Theo- 
rem 4.2) 
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y=ce" +e,¢* (5.6) 


In the special case 1, = —/,, the solution 5.6 can be rewritten as y = 
k, cosha,x+k, sinh /,x. 


Case 2. 1, =a+ib, a complex number. Since a, and ay in 5.1 and 
5.2 are assumed real, the roots of 5.2 must appear in conjugate pairs; 
thus, the other root is 2, = a — ib. Two linearly independent solutions 
are e4*")* and e")*, and the general complex solution is 


you en? + dee? (5.7) 


which is algebraically equivalent to 


y=ce™ cosbx+cye™ sin bx (5.8) 

Case 3. 4, =/,. Two linearly independent solutions are e** and 
xe*', and the general solution is 

y=ce*" +.c5xe"" (5.9) 


Warning: The above solutions are not valid if the differential equa- 
tion is not linear or does not have constant coefficients. Consider, for ex- 
ample, the equation y”—x’y =0. The roots of the characteristic equa- 
tion are A, =x and A, = —x, but the solution is not 


(x)x 


y=ceO* +c,e™ 


2 ) 
=cje" +ce~ 


Linear equations with variable coefficients are considered in Chapter 
Twelve. 


General Solution for nth-Order Equations 


The general solution of 5.3 is obtained directly from the roots of 5.4. If 
the roots 2,, A,,...,4,,are all real and distinct, the solution is 


y= ce +c,e" teetee (5.10) 
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If the roots are distinct, but some are com- 
plex, then the solution is again given by 5.10. 
As in the second-order equation, those terms 
involving complex exponentials can be com- 
bined to yield terms involving sines and co- 
sines. If A, is a root of multiplicity p [that is, if 
(A—4A,)? is a factor of the characteristic equa- 
tion, but (A — /,) P+! is not] then there will be 
p linearly independent solutions associated 
with A, given by et xe 70% xP 1e* These solutions are 
combined in the usual way with the solutions associated with the other 
roots to obtain the complete solution. 

In theory it is always possible to factor the characteristic equation, 
but in practice this can be extremely difficult, especially for differential 
equations of high order. In such cases, one must often use numerical tech- 
niques to approximate the solutions. See Chapter Fourteen. 


Solved Problems 
Solved Problem 5.1 Solve y”—y’—2y=0. 


The characteristic equation is 4? — 2 — 2 = 0, which can be factored 
into (A —2)(A + 1) = 0. Since the roots 1, = 2 and A, =—1 are real and dis- 
tinct, the solution is given by 5.6 as 


_ 2x -x 
y=cye + C7e 


Solved Problem 5.2 Solve y”—8y’+1l6y=0. 


The characteristic equation is 4” —8A+16=0 which can be fac- 
tored into (A — 4)? =0. The roots 2, =A, =4 are real and equal, so the gen- 
eral solution is given by 5.9 as 


y=ce +e)xe" 


Solved Problem 5.3 Solve y’” —6y”+2y’+36y=0. 


The characteristic equation A — 6A? + 22 + 36 = 0, has roots 2, = 
~2, A, =4+ivV2, and A, = 4—iV2.The solution is 
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y=ce 2 tdoetv2 4 yet 2): 


which can be rewritten, using Euler’s relations 


ib. 


e”* =cosbx+isinbx and e 


= cos bx —isin bx 


as 


y= Gee + dye* ei + dye*e 
y= Geo + dye (cos V2x+ isin ¥2x)+ d,e** (cos ¥2x —isin 2x) 
y=ce-* + (dy +d;)e*™ cos V2x + i(d, — d;)e** sin-/2x 


y=ce* +c,e" cos¥2x+c,e" sin V2x 


Note that this form of the solution corresponding to the complex roots 
can be easily formulated using Equation 5.8. 


Chapter 6 
SOLUTIONS 
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NONHOMOGENEOUS 
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The general solution to the linear differential equa- 


tion L(y) = ¢(x) is given by Theorem 4.4 as y= y, + 
y, where ys denotes one solution to the differential re \\ © 
equation and y, is the general solution to the asso- cA 

ciated homogeneous equation, L(y) = 0. Methods 
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for obtaining y, when the differential equation has constant coefficients 
are given in Chapter Five. In this chapter, we give methods for obtaining 
a particular solution y, once y, is known. 


The Method of Undetermined Coefficients 


Simple Form of the Method 


The method of undetermined coefficients is applicable only if @(x) and all 
of its derivatives can be written in terms of the same finite set of linearly 
independent functions, which we denote by {y,(x), y,(4),..., y,,(%)}. The 
method is initiated by assuming a particular solution of the form 


Vp (X) = Ay yy (%) + Ag yg (4) ++ + A, Vy (2X) 


where A,, A,,...,A,, denote arbitrary multiplicative constants. These arbi- 
trary constants are then evaluated by substituting the proposed solution 
into the given differential equation and equating the coefficients of like 
terms. 


Case 1. ¢(x) =p,,(x), an nth-degree polynomial in x. Assume a so- 
lution of the form 


yp =e (Ax +A, x be AL (6.1) 
where A, (j = 0,1,2,..., 2) is a constant to be determined. 
Case 2. ¢(x) = ke” where k and @ are known constants. Assume 
a solution of the form 

ae (6.2) 

where A is a constant to be determined. 
Case 3. @(x) =k, sin Bx + k, cos Bx where k,, k,, and # are known 
constants. Assume a solution of the form 


y= A sin Bx + B cos Bx (6.3) 


where A and B are constants to be determined. 
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Don’t Forget 


y= A sin Bx + B cos fx in its entirety is assumed 
for ¢(x) = k, sin Bx + k, cos Bx even when k, or k, 
is zero, because the derivatives of sines or cosines 
involve both sines and cosines. 


Generalizations 


If @(x) is the product of terms considered in Cases | through 3, take y, to 
be the product of the corresponding assumed solutions and algebraically 
combine arbitrary constants where possible. In particular, if @(x) =e p,,(x) 
is the product of a polynomial with an exponential, assume 


Vp HOM (A x” +A, ix | tet Ax t Ap) (6.4) 


where A. is as in Case 1. If, instead, @(x) = e** p_(x) sin Bx is the product 
of a polynomial, exponential, and sine term, or if @(x) = e™ p,,(x) cos Bx 
is the product of a polynomial, exponential, and cosine term, then assume 


Yp er i sin Bx(A,x” +A, 4x" + s+ Alyx + Ap) (6 5) 
+e™ cos Bx(B,x" +B, yx"! +---+B,x+By) 


where A; and B, (j = 0,1,2,..., ) are constants which still must be deter- 
mined. 
If @(x) is the sum (or difference) of terms already considered, then 
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we take y to be the sum (or difference) of the corresponding assumed so- 
lutions and algebraically combine arbitrary constants where possible. 


Modifications 

If any term of the assumed solution, disregarding multiplicative con- 
stants, is also a term of y, (the homogeneous solution), then the assumed 
solution must be modified by multiplying it by x”, where m is the small- 


est positive integer such that the product of x” with the assumed solution 
has no terms in common with y,. 


Limitations of the Method 


In general, if @(x) is not one of the types of functions considered above, 
or if the differential equation does not have constant coefficients, then the 
following method is preferred. 


Variation of Parameters 


Variation of parameters is another method for finding a particular solu- 
tion of the nth-order linear differential equation 


L(y) = 9(x) (6.6) 
once the solution of the associated homogeneous equation L(y) = 0 is 
known. Recall from Theorem 4.2 that if y,(), y,(),..., y,(x) are n linear- 
ly independent solutions of L(y) = 0, then the general solution of L(y) = 
Ois 


Va = AY (X) + Co yg (x) +++ +O, Yn (X) (6.7) 


The Method 
A particular solution of L(y) = ¢(x) has the form 


Vp = VY VoV2 Fo EVV y (6.8) 


where y, = y,(x) (i= 1,2,..., 2) is given in Equation 6.7 and v, (i= 1,2,..., 
n) is an unknown function of x which still must be determined. 
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To find v,, first solve the following linear equations simultaneously 
for vj: 


VY + VAY. t+ VAY, = 0 
VV +V3¥2 +--+ Vy, =O 


: (6.9) 
vil A ih teeta?) <0 


, -l , -l , -l 
Vi? VB VY? +--+ VED = O(a) 


Then integrate each v; to obtain v,, disregarding all constants of integra- 
tion. This is permissible because we are seeking only one particular so- 
lution. 


Example 6.1: For the special case n = 3, Equations 6.9 reduce to 


ViY, +V3Y2 +V393 = 0 
ViVi +V399 + V3¥3 = 0 (6.10) 
ViVi + V3y5 + V3y¥ = Ox) 


For the case n= 2, Equations 6.9 become 


Vi, + V3. =0 


org - , (6. 1 1 ) 
Vii + V2y2 = O(X) 
and for the case n= /, we obtain the single equation 
Vy, = (x) (6.12) 


Since y, (x), y5(%),..., y,,(x) aren linearly independent solutions of the same 
equation L(y) = 0, their Wronskian is not zero (Theorem 4.3). This means 
that the system 6.9 has a nonzero determinant and can be solved unique- 
ly for vj(x), V3 (X),...,¥), (x). 


Scope of the Method 


The method of variation of parameters can be applied to all linear differ- 
ential equations. It is therefore more powerful than the method of unde- 
termined coefficients, which is restricted to linear differential equations 
with constant coefficients and particular forms of ¢(x). Nonetheless, in 
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those cases where both methods are applicable, the method of undeter- 
mined coefficients is usually the more efficient and, hence, preferable. 

As a practical matter, the integration of v/(x) may be impossible to 
perform. In such an event other methods (in particular, numerical tech- 
niques) must be employed. 


Initial-Value Problems 


Initial-value problems are solved by applying the initial conditions to the 
general solution of the differential equation. It must be emphasized that 
the initial conditions are applied only to the general solution and not to 
the homogeneous solution y, that possesses all the arbitrary constants that 
must be evaluated. The one exception is when the general solution is the 
homogeneous solution; that is, when the differential equation under con- 
sideration is itself homogeneous. 


Solved Problems 


Solved Problem 6.1 Solve y”— y’—2y= 4x’. 


From Problem 5.1, y, = c,e* + c,e~. Here (x) = 4x", a second degree 
polynomial. Using Equation 6.1, we assume that 


y, =Agx? + Ax Ap (6.13) 


Thus, y, =2A,x+A, and y =2A,. Substituting these results into the 
differential equation, we have 


2A —(2Anx +A) —2(Agx” + Ayx +.Ag) = 4x7 
or, equivalently, 
(-2.Ay)x? +(-2A, — 2A, )x+(2A, — A, —2A)) =42x7 +(0)x +0 
Equating the coefficients of like powers of x, we obtain 


“GAp= de 2 9A2 S04 20) DAs oA 9A 0 
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Solving this system, we find that A, = —2, A, = 2, and A, = —3. Hence 
Equation 6.13 becomes 


Vy =-2x74+2x-3 
and the general solution is 


Y=, t+ Vp =Ce* tee * — 2x? +2x-3 


Solved Problem 6.2 Solve y”— y’—2y=sin2x. 

Again by Problem 5.1, y, =cje?* +c,e*. Here @(x) has the form dis- 
played in Case 3 with k, = 1, k, = 0, and B = 2. Using Equation 6.3, we 
assume that 


yp = Asin2x + Bcos2x (6.14) 


Thus, y,, = 2Acos2x—2Bsin2x and y; =—4Asin2x—4Bcos2x. Substi- 
tuting these results into the differential equation, we have 


(-4Asin2x —4Bcos2x)—(2Acos2x —2Bsin2.x) 
—2(Asin2x + Bcos2x) =sin2x 


or, equivalently, 
(-6A + 2B)sin2x +(-6B-—2A)cos2x = (1)sin2x +(0)cos2x 
Equating coefficients of like terms, we obtain 
-6A+2B=1 —2A —-6B=0 


Solving this system, we find that A = —3/20 and B = 1/20. Then from 
Equation 6.14, 


Yp =—Ssindx-+ = cos2x 


and the general solution is 
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y=ce +c,e° — er 2x+ ES cos 2x 
20 20 


Solved Problem 6.3 Solve y’”’ + y’=secx. 
This is a third-order equation with 
Vp = Cy + Cy COSX +3 SINX 
It follows from Equation 6.8 that 
Vp =V + V2 COS + V3 Sin x (6.15) 


Here y, = 1, y, =cos x, y, = sin x, and @(x) = sec x, so Equation 6.10 be- 
comes 
v1 (1) + v3 (cos x) + v3(sin x) = 0 
vj (0) + v3 (—sin x) + v4(cos x) =0 


v1 (0) + v5 (—cos x) + v3(—sin x) = sec x 


Solving this set of equations simultaneously, we obtain vj =secx, 
v, =—L, and v3, =—tanx. Thus, 


vy =| vide = Jsecxdx = Inlsecx + tan x | 


V> = [vide = [-Idx =-x 


sin x 


V3 = | vjde= J-tan xdx =—f ; 


dx = |n|cosx| 
Osx 


Substituting these values into Equation 6.15, we obtain 
V= In|secx+tanx|—xcosx+(sin x)In|cosx| 


The general solution is therefore 


Y=EYMtYp 


=cC, +C, cosx+c,sinx +In|secx + tan x|—xcos x + (sin x) In| cos x | 


Chapter 7 
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SECOND-ORDER 
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VL NNN8 


Spring Problems 


The simple spring system shown in Figure 7-1 consists of a mass m at- 
tached to the lower end of a spring that is itself suspended vertically from 
a mounting. The system is in its equilibrium position when it is at rest. 
The mass is set in motion by one or more of the following means: dis- 
placing the mass from its equilibrium position, providing it with an ini- 
tial velocity, or subjecting it to an external force F(#). 
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Equilibrium position Initial position at t= 0 


Y 


Y 


x=0 


F(t) 


Positive x-direction 


Figure 7-1 


Hooke’s law: The restoring force F of a spring is equal and opposite to 
the forces applied to the spring and is proportional to the extension (con- 
traction) | of the spring as a result of the applied force; that is, F = —kl, 
where k denotes the constant of proportionality, generally called the 
spring constant. 


Example 7.1. A steel ball weighing 128 Ib is suspended from a spring, 
whereupon the spring is stretched 2 ft from its natural length. The applied 
force responsible for the 2-ft displacement is the weight of the ball, 128 Ib. 
Thus, F = —128 lb. Hooke’s law then gives —128 = —k(2), or k = 64 lb/ft. 


For convenience, we choose the downward direction as the positive 
direction and take the origin to be the center of gravity of the mass in the 
equilibrium position. We assume that the mass of the spring is negligible 
and can be neglected and that air resistance, when present, is proportion- 
al to the velocity of the mass. Thus, at any time ¢, there are three forces 
acting on the system: (1) F(t), measured in the positive direction; (2) a 
restoring force given by Hooke’s law as F’, = —kx, k > 0; and (3) a force 
due to air resistance given by F’, = —ax, a > 0, where a is the constant of 
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proportionality. Note that the restoring force F always acts in a direction 
that will tend to return the system to the equilibrium position: if the mass 
is below the equilibrium position, then x is positive and —kx is negative; 
whereas if the mass is above the equilibrium position, then x is negative 
and —kx is positive. Also note that because a > 0 the force F’, due to air 
resistance acts in the opposite direction of the velocity and thus tends to 
retard, or damp, the motion of the mass. 

It now follows from Newton’s second law (see Chapter Three) that 
mx = —kx — ax + F(t), or 


fp Ot ge (7.1) 
mom m 
If the system starts at ¢= 0 with an initial velocity v, and from an initial 


position x), we also have the initial conditions 


The force of gravity does not explicitly appear in 7.1, but it is pres- 
ent nonetheless. We automatically compensated for this force by mea- 
suring distance from the equilibrium position of the spring. If one wish- 
es to exhibit gravity explicitly, then distance must be measured from the 
bottom end of the natural length of the spring. That is, the motion of a vi- 
brating spring can be given by 

el Bs F(t) 
X+—x+—x=g+— 
m om 
if the origin, x = 0, is the terminal point of the unstretched spring before 
the mass m is attached. 


Electrical Circuit Problems 


The simple electrical circuit shown in Figure 7-2 consists of a resistor R 
in ohms; a capacitor C in farads; an inductor L in henries; and an elec- 
tromotive force (emf) E(f) in volts, usually a battery or a generator, all 
connected in series. The current / flowing through the circuit is measured 
in amperes and the charge g on the capacitor is measured in coulombs. 
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+ 
E(t) Cc 


Figure 7-2 


Kirchhoff’s loop law: The algebraic sum of the voltage drops in a sim- 
ple closed electric circuit is zero. 


It is known that the voltage drops across a resistor, a capacitor, and 
an inductor are respectively RI, (1/C)q, and L(d//dt) where q is the charge 
on the capacitor. The voltage drop across an emf is —E(t). Thus, from 
Kirchhoff’s loop law, we have 


dl 1 
RI+L—+—q-E(t)=0 7.3 
anc! (t) (7.3) 
The relationship between g and J is 
2 
7a 4 ea (7.4) 
dt dt dt 
Substituting these values into Equation 7.3, we obtain 
d’q Rdq 1 1 
+ + =— E(t 75 
eae be ae oP 
The initial conditions for g are 
d 
q(9) = do <1 =1(0)=1 (7.6) 
dt |,-0 


To obtain a differential equation for the current, we differentiate 
Equation 7.3 with respect to ¢ and then substitute Equation 7.4 directly 
into the resulting equation to obtain 

d’I Rdl 1 1 dE(t) 
>t +—J= 
dt’ Ldt LC L at 


(7.7) 
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The first initial condition is 1(0) = J,. The second initial condition is ob- 
tained from Equation 7.3 by solving for d//dt and then setting t= 0. Thus, 
dl 1 R 1 
=—E(0 I 7.8 
Hills Fe L° ic 1° (7.8) 
An expression for the current can be gotten either by solving Equation 
7.7 directly or by solving Equation 7.5 for the charge and then differ- 
entiating that expression. 


Buoyancy Problems 


Consider a body of mass m submerged either partially or totally in a liq- 
uid of weight density p. Such a body experiences two forces, a downward 
force due to gravity and a counter force governed by: 


Archimedes’ principle: A body in liquid experiences a buoyant upward 
force equal to the weight of the liquid displaced by that body. 


Positive x-direction 


Equilibrium state Equilibrium position 


i 
x(t) 
: ; a ios it is ee 


Water line 


Figure 7-3 


Equilibrium occurs when the buoyant force of the displaced liquid 
equals the force of gravity on the body. Figure 7-3 depicts the situation 
for a cylinder of radius r and height H where h units of cylinder height 
are submerged at equilibrium. At equilibrium, the volume of water dis- 
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placed by the cylinder is zr7h, which provides a buoyant force of zr7hp 
that must equal the weight of the cylinder mg. Thus, 


ar-ho = mg (7.9) 
Motion will occur when the cylinder is displaced from its equilibri- 
um position. We arbitrarily take the upward direction to be the positive 
x-direction. If the cylinder is raised out of the water by x(f) units, as shown 
in Figure 7-3, then it is no longer in equilibrium. The downward or neg- 
ative force on such a body remains mg but the buoyant or positive force 
is reduced to zr*[h — x(t)]p. It now follows from Newton’s second law 
that 
mx = mr? [h— x(t)]o —mg 
Substituting 7.9 into this last equation, we can simplify it to 


mx = —mr’x(t)p 


or 


2 
fe Pye (7.10) 


Classifying Solutions 


Vibrating springs, simple electrical circuits, and floating bodies are all 
governed by second-order linear differential equations with constant co- 
efficients of the form 


X+axt+ayx = f(t) (7.11) 


For vibrating spring problems defined by Equation 7.1, a, = a/m, a) = 
k/m, and f(t) = F()/m. For buoyancy problems defined by Equation 7.10, 
a, =0, a, = a1’p/m, and f(t) = 0. For electrical circuit problems, the in- 
dependent variable x is replaced either by g in Equation 7.5 or J in Equa- 
tion 7.7. 

The motion or current in all of these systems is classified as free and 


CHAPTER 7: Second-Order Linear Differential Equations 53 


undamped when f(t) = 0 and a, = 0. It is classi- 
fied as free and damped when f(t) is identically 
zero but a, is not zero. For damped motion, 
there are three separate cases to consider, ac- 
cording as the roots of the associated character- 
istic equation (see Chapter Five) are (1) real and 
distinct, (2) equal, or (3) complex conjugate. 
These cases are respectively classified as (1) 
overdamped, (2) critically damped, and (3) os- 
cillatory damped (or, in electrical problems, underdamped ). If f(t) is not 
identically zero, the motion or current is classified as forced. 

A motion or current is transient if it “dies out” (that is, goes to zero) 
as t > ce. A steady-state motion or current is one that is not transient and 
does not become unbounded. Free damped systems always yield transient 
motions, while forced damped systems (assuming the external force to be 
sinusoidal) yield both transient and steady-state motions. 

Free undamped motion defined by Equation 7.11 with a, = 0 and 
f() = 0 always has solutions of the form 


x(t) = c, cos@t + cy sin@t (7.12) 


which defines simple harmonic motion. Here c,, c,, and @ are constants 


with @ often referred to as circular frequency. The natural frequency f is 


and it represents the number of complete oscillations per time unit un- 
dertaken by the solution. The period of the system of the time required to 
complete one oscillation is 


Tas 
z 


Equation 7.12 has the alternate form 
x(t) = (-1)* Acos(at — ¢) (7.13) 


where the amplitude A= Kien +05 , the phase angle ¢ = arctan(c,/c,), 
and k is zero when c, is positive and unity when c, is negative. 
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Solved Problems 


Solved Problem 7.1 A 10-kg mass is attached to a spring, stretching it 
0.7 m from its natural length. The mass is started in motion from the 
equilibrium position with an initial velocity of | m/sec in the upward di- 
rection. Find the subsequent motion, if the force due to air resistance is 
—90x N. 


Taking g = 9.8m/sec”, we have w = mg = 98 N and k= w/l= 140 N/m. 
Furthermore, a = 90 and F(t) =0 (there is no external force). Equation 7.1 
becomes 


X¥+9x+14x=0 (7.14) 


The roots of the associated characteristic equation are A, = —2 and A, = 
—7, which are real and distinct; hence this problem is an example of over- 
damped motion. The solution of 7.14 is 


x= eer + cet 


The initial conditions are x(0) = O (the mass starts at the equilibrium 
position) and x(0) = —1 (the initial velocity is in the negative direction). 


Applying these conditions, we find that c, =—c, =-4t , so that x = 


t(e™ ey, Note that x > 0 as t  o9; thus, the motion is transient. 


Solved Problem 7.2 Determine whether a cylinder of radius 4 in, height 
10 in, and weight 15 Ib can float in a deep pool of water of weight densi- 
ty 62.5 Ib/ft?. 


Let h denote the length (in feet) of the submerged portion of the 
cylinder at equilibrium. With r= + ft, it follows from Equation 7.9 that 


mg _ 15 
mrp  m(4)°62.5 


= 0.688ft = 8.25in 


Thus, the cylinder will float with 10 — 8.25 = 1.75 in of length above the 
water at equilibrium. 


Chapter 8 
LAPLACE 


TRANSFORMS AND 
INVERSE LAPLACE 
TRANSFORMS 


In THIS CHAPTER: 


Y Definition of the Laplace Transform 
Y Properties of Laplace Transforms 
Y Definition of the Inverse Laplace 
Transform 

Manipulating Denominators 
Manipulating Numerators 
Convolutions 

Unit Step Function 

Translations 

Solved Problems 


CNV s 


55 


Copyright 2003 by The McGraw-Hill Companies, Inc. Click Here for Terms of Use. 


56 DIFFERENTIAL EQUATIONS 
Definition of the Laplace Transform 
Let f(x) be defined for 0 < x < © and let s denote an arbitrary real vari- 


able. The Laplace transform of f(x), designated by either & { f(x)} or F(s), 
is 


SL f(x)} = F(s) = fe" f(a)dx (8.1) 
0 


for all values of s for which the improper integral converges. Conver- 
gence occurs when the limit 


R 
lim fer f(x)dx (8.2) 
Roo 0 


exists. If this limit does not exist, the improper integral diverges and f(x) 
has no Laplace transform. When evaluating the integral in Equation 8.1, 
the variable s is treated as a constant because the integration is with re- 
spect to x. 

The Laplace transforms for a number of elementary functions are 
given in Appendix A. 


Properties of Laplace Transforms 
Property 8.1 (Linearity). If £{ f(x)} = F(s) and L{g(x)} = G(s), then 
for any two constants c, and c, 

Lie, f@) + c,g@)} =c, L{f@)} + co L{gQ)} (8.3) 

=c,F(s) + c,G(s) 
Property 8.2. If £{ f(x)} = F(s), then for any constant a 
Lf{e“f(x)} = F(s — a) (8.4) 

Property 8.3. If £{ f(x)} = F(s), then for any positive integer n 


d” 
F 
as [F(s)] (8.5) 


L{x"f(x)} = (-D" 
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Property 8.4. If £{ f(x)} = F(s) and if _ fe) exists, then 
x00 X 
x>0 
i S 
{* fc} = | F(nat (8.6) 
x Ss 
Property 8.5. If £{ f(x)} = F(s), then 
t 1 
aj so} = — F(s) (8.7) 
s 


Property 8.6. If f(x) is periodic with period a, that is, f(x + @) =f(x), then 


fe“ s@)ax 
Lf (x)}= arr anl (8.8) 


Functions of Other Independent Variables 


For consistency only, the definition of the Laplace transform and its prop- 
erties, Equations 8.1 through 8.8, are presented for functions of x. They 
are equally applicable for functions of any independent variable and are 
generated by replacing the variable x in the above equations by any vari- 
able of interest. In particular, the counter part of Equation 8.1 for the 
Laplace transform of a function of f is 


SA f(t)} = F(s)= fe fat 
0 


Definition of the Inverse Laplace Transform 


An inverse Laplace transform of F(s) designated by £-!{ F(s)}, is an- 
other function f(x) having the property that £{ f(x)} = F(s). 

The simplest technique for identifying inverse Laplace transforms is 
to recognize them, either from memory or from a table such as in the Ap- 
pendix. If F(s) is not in a recognizable form, then occasionally it can be 
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transformed into such a form by algebraic manipulation. Observe from 
the Appendix that almost all Laplace transforms are quotients. The rec- 
ommended procedure is to first convert the denominator to a form that 
appears in the Appendix and then the numerator. 


Manipulating Denominators 


The method of completing the square converts a quadratic polynomial 
into the sum of squares, a form that appears in many of the denominators 
in the Appendix. In particular, for the quadratic 


b 
as’ +bs+c =o 5? +25 +C 


=a(s+k) +h? 


where k = b/2a and h=4\c—(b* / 4a). 


The method of partial fractions transforms a function of the form 
a(s)/b(s), where both a(s) and b(s) are polynomials in s, into the sum of 
other fractions such that the denominator of each new fraction is either a 
first-degree or a quadratic polynomial raised to some power. The method 
requires only that the degree of a(s) be less than the degree of b(s) (if this 
is not the case, first perform long division, and consider the remainder 
term) and b(s) be factored into the product of distinct linear and quadrat- 
ic polynomials raised to various powers. 

The method is carried out as follows. To each factor of b(s) of the 
form (s — a)", assign a sum of m fractions, of the form 


A A, An 
+ zt t m 
s-a_ (s—a) (s—a) 


To each factor of b(s) of the form (s* + bs + c)?, assign a sum of p frac- 
tions, of the form 


Bys+C, mn Bys+C, M Bs+C, 
S+bstc (s? +bs+c)? (s? +bs+c)? 
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Here A,, B, and C, (= 1,2,..., m; j, k= 1,2,..., p) are constants which still 
must be determined. 

Set the original fraction a(s)/b(s) equal to the sum of the new frac- 
tions just constructed. Clear the resulting equation of fractions and then 
equate coefficients of like powers of s, thereby obtaining a set of simul- 
taneous linear equations in the unknown constants A.,, B, and C,,. Final- 
ly, solve these equations for A,, B, and C,. 


Manipulating Numerators 
A factor s — a in the numerator may be written in terms of the factor s — 
b, where both a and b are constants, through the identity s—a=(s— b)+ 


(b—a). The multiplicative constant a in the numerator may be written ex- 
plicitly in terms of the multiplicative constant b through the identity 


a 
G= (0) 


Both identities generate recognizable inverse Laplace transforms when 
they are combined with: 


Property 8.7 (Linearity). If the inverse Laplace transforms of two func- 
tions F(s) and G(s) exist, then for any constants c, and c,, 


Lc, F(s)+cyG(s)}=c, £7 {F(s)} +c, £7" {G(s)} 


Convolutions 


The convolution of two functions f(x) and g(x) is 


f(x)* g(a) = J fg(e—Nar 
0 (8.9) 


Theorem 8.1. f(x) * g(x) = g(x) * f(x). 


Theorem 8.2. (Convolution Theorem). If £{f(x)} = F(s) and £{g(x)} 
= G(s), then L{ f(x) * g(x)} = L{f)} L{g(a)} = F(s)G(s) 


60 DIFFERENTIAL EQUATIONS 


You Need to Know S 


The inverse Laplace transform of a product is 
computed using a convolution. 


LF (s)G(s)} =f(x) *9(x) = g(x) *f(x) (8.10) 


If one of the two convolutions in Equation 8.10 is simpler to calculate, 
then that convolution is chosen when determining the inverse Laplace 
transform of a product. 


Unit Step Function 


The unit step function u(x) is defined as 


7 0 x<0 
Us 1 x20 


As an immediate consequence of the definition, we have for any num- 
ber c, 


0 x<c 
WO) 1 x2c 


The graph of u(x — c) is given in Figure 8-1. 


u(x—c) 


Figure 8-1 
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Theorem 8.3. L{u(x—c)}= : er. 
s 


Translations 
Given a function f(x) defined for x = 0, the function 


0 xXx<C 
f(x-c) x2c 


Wx-ofs-0=| 


represents a shift, or translation, of the function f(x) by c units in the pos- 
itive x-direction. For example, if f(x) is given graphically by Figure 8-2, 
then u(x — c)f(x — c) is given graphically by Figure 8-3. 


fix) 


Figure 8-2 


udx-c)fixteoc) 


Figure 8-3 
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Theorem 8.4. If F(s) = L£{ f(x}, then 
L{u(x—c)f(x-c)} =e % F(s) 


Conversely, 


x<c 
f(x-c) x2c 


Ee" FG) = Ws-of-9 =| 


Solved Problems 
Solved Problem 8.1 Find £{e”}. 


Using Equation 8.1, we obtain 


2 R 
F(s)= feed = jim feordx 


0 0 
=R 
fF elaox P Pele 4 
= lim = lim | ————— 
Roo a—s Reo a-s 
x=0 
= (for s > a) 
s-—a 


Note that when s < a, the improper integral diverges. (See also entry 7 in 
the Appendix.) 


Solved Problem 8.2 Find £{xe**}. 
This problem can be done three ways. 


(a) Using entry 14 of the Appendix with n = 2 and a = 4, we have 
directly that 
1 


yf 4x = 
{xe""} (aay 


(b) Set f(x) =x. Using Property 8.2 with a=4 and entry 2 of the Ap- 
pendix, we have 
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1 
F(s)=L£{f(x)}=L{x} ==> 


Ss 


and 
4x = = = —_| 
L{e™x}= F(s—4) GD: 


(c) Set f(x) = e*. Using Property 8.3 with n = 1 and the results of 
Problem 8.1, or alternatively, entry 7 of the Appendix with a= 4 
we find that 


1 
s—4 


F(s)= L{ f(x) } = Lle*} = 


and 


Axio te he) af RoW ¢ 2 
ae eee (4) (s—4)" 
s+3 


Solved Problem 8.3 Use partial fractions to decompose ——————-. 
(s—2)(s+1) 


To the linear factors s — 2 and s + 1, we associate respectively the 
fractions A/(s — 2) and B/(s + 1). We set 


s+3 A B 


a + — 
(s—2)(st+1l) s—-2 stl] 


and, upon clearing fractions, obtain 


s+3=A(s+1)+ B(s—2) (8.11) 


To find A and B, we substitute s =—1 and s = 2 into 8.11, we immediate- 
ly obtain A = 5/3 and B = —2/3. Thus, 


s+3 _ 5/3 2/3 


(s—2)(s+1) s—2 stl 


2 3 
a atcepeal yt _ E Ys 22 | 
Solved Problem 8.4 Find fe “D+ 
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No function of this form appears in the Appendix. Using the results 
of Problem 8.3 and Property 8.7, we obtain 


wife Sarl Uh Ba) 
(s—2)(s+){ 3 §=2) 3 s+l 


=> o2x_2 -—x 
3 3 
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Laplace Transforms of Derivatives 


Denote £{y(x)} by Y(s). Then under very broad conditions, the Laplace 
transform of the nth-derivative (n = 1,2,3,...) of y(x) is 


n Pe! _ gal _ gt-2 ys i 
gy d u =s Y(s)-s *) Ss yy (9.1) 
dx Say Oy 0) 
If the initial conditions on y(x) at x = 0 are given by 
y(0) = Cp, yO) = C15--¥" (0) = C44 (9.2) 
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then (9.1) can be rewritten as 


d” n n- n- 
“| o} =s"Y(s)— cys" — C8" 7 Cy 98 — Gy, (9.3) 
For the special cases of n = 1 and n = 2, Equation 9.3 simplifies to 


L{y’(x)} = sY(s)— co (9.4) 


Lly"(x)} = 8°¥(s)—egs— cy (9.5) 


6 Note! 


Laplace transforms convert differential equations 
into algebraic equations. 


Solutions of Linear Differential Equations 
with Constant Coefficients 


Laplace transforms are used to solve initial-value problems given by the 
nth-order linear differential equation with constant coefficients 
n-1 
eae ane AD each © eb ye) (9.6) 
dx dx 
together with the initial conditions specified in Equation 9.2. First, take 
the Laplace transform of both sides of Equation 9.6, thereby obtaining an 
algebraic equation for Y(s). Then solve for Y(s) algebraically, and final- 
ly take inverse Laplace transforms to obtain y(x) = £-!{ ¥(s)}. 

Unlike previous methods, where first the 
differential equation is solved and then the 
initial conditions are applied to evaluate 
the arbitrary constants, the Laplace trans- 
form method solves the entire initial-value 
problem in one step. There are two excep- 
tions: when no initial conditions are speci- 
fied and when the initial conditions are not 
at x = 0. In these situations, cy through c,,_, 
in Equations 9.2 and 9.3 remain arbitrary 
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and the solution to differential equation 9.6 is found in terms of these con- 
stants. They are then evaluated separately when appropriate subsidiary 
conditions are provided. 


Solutions of Linear Systems 


Laplace transforms are useful for solving systems of linear differential 
equations; that is, sets of two or more differential equations with an equal 
number of unknown functions. If all of the coefficients are constants, then 
the method of solution is a straightforward generalization of the one de- 
scribed above. Laplace transforms are taken of each differential equation 
in the system; the transforms of the unknown functions are determined 
algebraically from the resulting set of simultaneous equations; inverse 
transforms for the unknown functions are calculated with the help of the 
Appendix. 


Solved Problems 
Solved Problem 9.1 Solve y’—5y=e>*; y(0)=0. 
Taking the Laplace transform of both sides of this differential equation 


and using Property 8.1, we find that L{y’}-5 ¥F {y}= Lf{e>*}. Then, us- 
ing the Appendix and Equation 9.4 with c, = 0, we obtain 


1 : 1 
[sY(s)—0]—5Y(s) 5-5 from which Y(s)= Gos? 


Finally, taking the inverse Laplace transform of Y(s), we obtain 


ie _ 1 
yxy=LUY(S)}=L |} nes 


(see Appendix, entry 14). 


Solved Problem 9.2 Solve the system 
y"+z+y=0 
zt+y=0; 
y(0)=0, y'(0)=0, z2(0)=1 
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Denote L£{y(x)} and L{z(x)} by Y(s) and Z(s) respectively. Then, 
taking the Laplace transforms of both differential equations, we obtain 


[s°¥(s)—(0)s—(0)]+ Z(s) + ¥(s)=0 
[sZ(s)—1]+[sY(s)—0]=0 


or 


(s? +D)Y(s)+ Z(s)=0 


Y(s)+Z(s) = Es 
s 


Solving this last system for Y(s) and Z(s), we find that 


1 11 
Y(s)=-— Z(s)=-+5 
S Ss AY 


Thus, taking inverse transforms, we conclude that 


1 1 
Wa)a-5e ax)=1+5x° 
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Matrices and Vectors 


A matrix (designated by an uppercase boldface letter) is a rectangular ar- 
ray of elements arranged in horizontal rows and vertical columns. In this 
book, the elements of matrices will always be numbers or functions of 
the variable ¢. If all the elements are numbers, then the matrix is called a 
constant matrix. 


Example 10.1. 


1 2/|1 e 2 
R ,and}]1 t* cost 
3 4}|/¢ -1 1 
are all matrices. In particular, the first matrix is a constant matrix, where- 


as the last two are not. 
A general matrix A having p rows and n columns is given by 


a, Ap Ay 

ay, 42 Ayn, 
A= [a,| 

Ay Ayo Oo 


where a.. represents that element appearing in the ith row and jth column. 
a 


he Note! 


Amatrix is square if it has the same number of rows 
and columns. 


A vector (designated by a lowercase boldface letter) is a matrix having 
only one column or one row. (The third matrix given in Example 10.1 is 
a vector.) 
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Matrix Addition 


The sum A + B of two matrices A = la,/] and B = [b;,] having the same 
number of rows and the same number of columns is the matrix obtained 
by adding the corresponding elements of A and B. That is, 


A+B=[a;]+[bj;]=[aj + 5;] 


Matrix addition is both associative and commutative. Thus, A+(B+C) 
=(A+B)+C and A+B=B+A. 


Scalar and Matrix Multiplication 


If A is either a number or a function of t, then AA (or, equivalently, AA) is 
defined to be the matrix obtained by multiplying every element of A by 
A. That is, 


AA = Aa; |= [Aaj] 


Let A= [a,,] and B= [b;,] be two matrices such that A has r rows and n 
columns and B has n rows and p columns. Then the product AB is de- 
fined to be the matrix C = [c;] given by 


= ab. Gal.) 
k=1 


The element c,,is obtained by multiplying the elements of the ith row of 
A with the corresponding elements of the jth column of B and summing 
the results. 

Matrix multiplication is associative and distributes over addition; 
in general, however, it is not commutative. Thus, 


A(BC) = (AB)C, A(B + C) = AB+ AC, and (B+ C)A=BA+ CA. 


Important! 


In general, AB + BA. 
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Powers of a Square Matrix 


If n is a positive integer and A is a square matrix, then 


A" =AA---A 
—— 


n_ times 


In particular, A? = AA and A? = AAA. By definition, A° = I, where 


1 0 0 0 0 

0 1 0 0 0 

pe 0 0 1 0 : 
0 

L 0 1| 


is called an identity matrix. For any square matrix A and identity matrix 
I of the same size 


AI=IA=A 


Differentiation and Integration of Matrices 


The derivative of A= la;,] is the matrix obtained by differentiating each 


element of A; that is, 
a8) 5 
dt dt 


Similarly, the integral of A, either definite or indefinite, is obtained by in- 
tegrating each element of A. Thus, 


a 


b 
JAar= 


J ot 


ja and [ Adr= 
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The Characteristic Equation of a Matrix 


The characteristic equation of a square matrix A is the polynomial equa- 
tion in A given by 


det(A — AI) = 0 (10.1) 
where det() stands for “the determinant of.” Those values of 4 which sat- 
isfy 10.1, that is, the roots of 10.1, are the eigenvalues of A, a k-fold root 


being called an eigenvalue of multiplicity k. 


Theorem 10.1. (Cayley-Hamilton theorem). Any square matrix satis- 
fies its own characteristic equation. That is, if 


det(A — AI) =b,A" +b,_)A" | +--+ BA? +b,A + by 


then b,A" +b, A") +--+b,A? +b At+boI =0. 


n-1 
Definition of the Matrix Exponential e4! 
For a square matrix A, 

Melts “Ar MP Dr —AN (10.2) 


The infinite series 10.2 converges for every A He t, so that e“? is defined 
for all square matrices. 


Computation of the Matrix Exponential e4¢ 


For actually computing the elements of e“’, 10.2 is not generally useful. 
However, it follows (with some effort) from Theorem 10.1, applied to the 
matrix Af, that the infinite series can be reduced to a polynomial in ¢. 
Thus: 


Theorem 10.2. If A is a matrix having n rows and n columns, then 


ea, AT +a, Ate ++ a,Att+agl (10.3) 
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where O15, O,,...,0! are functions of t which must be determined for 


each A. ite 
Example 10.2. When A has two rows and two columns, then n = 2 and 
e“ =a,At+aol (10.4) 
When A has three rows and three columns, then n = 3 and 
e“' =a,A°? +a,At+aol (10.5) 
Theorem 10.3. Let A be as in Theorem 10.2, and define 
(A) =A" | +d” 2 +++ aA? +a,A4+ a (10.6) 
Then if A, is an eigenvalue of At, 
e =r(A,) (10.7) 


Furthermore, if 2, is an eigenvalue of multiplicity k, k > 1, then the fol- 
lowing equations are also valid: 


et = oa 
de aan 
z 
efi —zr(A) 
a es (10.8) 
k-1 
e parr) 
A=A, 


Note that Theorem 10.3 involves the eigenvalues of Af; these are t 
times the eigenvalues of A. When computing the various derivatives in 
10.8, one first calculates the appropriate derivatives of the expression 
10.6 with respect to A, and then substitutes 2 = 1,. The reverse procedure 
of first substituting 2 = 2, (a function of f) into 10.6, and then calculating 
the derivatives with respect to f can give erroneous results. 
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Example 10.3. Let A have four rows and four columns and let 2 = 5t and 
A = 2t be eigenvalues of At of multiplicities three and one, respectively. 
Then n= 4 and 
r(A)= a3a3 + yA” +ajA+a 
r'(A) = 304A" +2a,A +a, 
r’(A) = 6031 + 2a, 


Since A = 5t is an eigenvalue of multiplicity three, it follows that e>’ = 
r(5t), e? = r' (52) and e* = r'(5t). Thus, 

e” =a, (5t)? +a (St) +a, (5t)+ a9 
e” =3a,(5t) + 2a,(5t) +a, 


e =6a;(5t) + 2a, 
Also, since A = 2t is an eigenvalue of multiplicity one, it follows that e’ 
= r(2ft), or 


e” =a3(2t) +a5(2t) +a, (21) +a 


Notice that we now have four equations in four unknown a’s. 


Method of computation: For each eigenvalue A, of At, apply Theorem 
10.3 to obtain a set of linear equations. When this is done for each eigen- 
value, the set of all equations so obtained can be solved for a,q,,,..., 
@,,_,. These values are then substituted into Equation 10.3 which, in turn, 
is used to compute e4*. 


Solved Problems 
Solved Problem 10.1 Find AB and BA for 
12 3 7 O 
A= R B= E 
laser Bla al 


Since A has three columns and B has two rows, the product AB is not de- 
fined. But 


76 DIFFERENTIAL EQUATIONS 


[7 O]1 2 3 
BA 

js ils 5 ol 
_| 7+ OA) 7(2) + (0)(5) heed 
L8d) +(-D(4) 8(2)+(-DG) 83)+(-D) 


[7 14 21 
(4 11 18 


1 3 
Solved Problem 10.2 Find the eigenvalues of A = ; } 


We have 
1 3 1 O 
A-Al= +(-A) 
4 2 0 1 
“4 1 3 n —A 0 ie 1-A 3 
See Oe aa a een 
1-A 3 
det(A — AD = ce 
Hence, 4 2-a 


=(1—A)(2—A)—(3)(4) = A? —34-10 


The characteristic equation of A is A* — 3A —10=0, which can be factored 
into (A — 5)(A + 2) = 0. The roots of the equation are 1, = 5 and 2, = —2, 
which are the eigenvalues of A. 


1 1 
Solved Problem 10.3 Find e“! for A = i i! 


Here n = 2. From Equation 10.4, 


(10.9) 


att+aAg at 
gS a, At+ aol = 


9ayt ayt + A 
and from Equation 10.6, r(A) = @,A + a. The eigenvalues of Af are 1, = 


4t and A, =—21, which are both of multiplicity one. Substituting these val- 
ues successively into Equation 10.7, we obtain the two equations 


e' = dta, +a 


e! = —2tar, + ag 
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Solving these equations for @, and a, we find that 


= =(e" ag) and a= ste + 2e7') 


Substituting these values into 10.9 and simplifying, we have 


ee oi Fe 43e72! ett — et | 
ett = 


6| 9e4 902! 3e4 4. 3e72! 


0 1 
Solved Problem 10.4 Find e4‘ for A= |: ‘| 


Since n = 2, it follows from Equations 10.4 and 10.6 that 
a ant 
eM =a Attagl=| f ‘ (10.10) 
8ayt —2ayt + Ao 
and r(A) = a,A + a. The eigenvalues of At are 2, = 2tand A, =—41, which 


are both of multiplicity one. Substituting these values successively into 
Equation 10.7, we obtain 


e* =a, (2t) + a 


e' =a,(-4t) +a 


Solving these equations for @, and a, we find that 


= ae — ay and a= (oem + ey 


Substituting these values into 10.10 and simplifying, we have 


Ar_ 1 Ae Ge ag oee 


6:82" =8e * Bert 44e™ 
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Reduction of Linear Differential Equations 
to a First-Order System 


Reduction of One Equation 
Every initial-value problem of the form 
d"x a ly 
b, oo ae Os aa “tbh (Hx+bo(x= g(t); CLD 
Mi jN=y.° AGS t pack” yates (11.2) 
with b(t) # 0, can be reduced to the first-order matrix system 


x(t) = A(t)x(t) +f (1) 


X(t) ) =e 


(11.3) 


where A(?), f(¢), ¢, and the initial time 7, are known. The method of re- 
duction is as follows. 


Step 1. Rewrite 11.1 so that d’x/dt" appears by itself. Thus, 


Mx 4 tt 
dt” = ay-| dt n-1 (11.4) 
+e +ay(t)x+ag(t)x+ f(t) 


where aft) = —b{tlb,(t) (j=0,1,...,— 1) and fd) = g(1)/b(). 


Step 2. Define 1 new variables (the same number as the order of the orig- 
inal differential equation), x,(¢), x,(2),..., x,(0), by the equations 


BO=eO, BOs aa 
P “ ae és (11.5) 
x(t x(t 
X3(t)= ae pee X, 1) = dt 


These new variables are interrelated by the equations 
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X(t) = X(t) 
Xy (t)= X3 (t) 
X3(t) = x4(t) (11.6) 


X,-1(t) = X, (0) 


Step 3. Express d*/dt in terms of the new variables. Proceed by first dif- 
ferentiating the last equation of 11.5 to obtain 


d Pa _ d"x(t) 


An(f) dt| ar"! dt” 


Then, from Equations 11.4 and 11.5, 


n-1 
x,@) =a,4 Ce +. +a, (t)x(t)+ dg (t)x()+ f@) 


= Ay (EX, (1) ++. (t)X9(t) + ag (x (+ fF) 


For convenience, we rewrite this last equation so that x (4) appears before 
x,(t), etc. Thus, 


Xp (1) = do (tx (1) + a (x2 (t)h+--- +a, (Ox, + f(t) C17) 


Step 4. Equations 11.6 and 11.7 are a system of first-order linear differ- 
ential equations in x, (1), x,(0),..., x,(0). This system is equivalent to the 
single matrix equation x(f) = A(f)x(t) + f(® if we define 


x, (t) 


X(t) 


x(t) = (11.8) 


x, (t) 


f()=| : (11.9) 
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0 1 0 O  -: 0 
0 0 0 0 
0 0 0 1 tee 0 
A(t)= . . ; . . (11.10) 
0 0 0 O  -: 1 
| do (t) a(t) ay (t) as (t) a, (t) | 


Step 5. Define 


Then the initial conditions (11.2) can be given by the matrix (vector) 
equation x(f,) = c. This last equation is an immediate consequence of 
Equations 11.8, 11.5, and 11.2, since 


X; (to) X(t) co 


2052 ay _ Xf) = a 


Ml 
oO 


como ie ke ation 67-9 4 pone 
Observe that if no initial conditions are prescribed, Steps | through 4 by 
themselves reduce any linear differential equation 11.1 to the matrix 
equation X(t) = A(x(f) + f(A). 


Reduction of a System 


Aset of linear differential equations with initial conditions also can be re- 
duced to system 11.3. The procedure is nearly identical to the method for 
reducing a single equation to matrix form; only Step 2 changes. With a 
system of equations, Step 2 is generalized so that new variables are de- 
fined for each of the unknown functions in the set. 
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Solution of the Initial-Value Problem 


By the procedure described above, any initial-value problem in which the 
differential equations are all linear with constant coefficients, can be re- 
duced to the matrix system 


X(t) =Ax(1)+f();  x(ty) =e (11.11) 


where A is the matrix of constants. The solution to Equation 11.11 is 


t 
x(t) = Mes eM [eM E(s)ds (11.12) 


to 


or equivalently 


t 
x(t) = ee | eA F(s)ds (11.13) 


to 


In particular, if the initial-value problem is homogeneous [i.e., f(t) = 0], 
then both equations 11.12 and 11.13 reduce to 


x(t) = eC We (11.14) 


A(t-to) 3-As A(t— 
> 


In the above solutions, the matrices e ,e and e4¢-») are easily 
computed from e“’ by replacing the variable ¢ by t — f,, —s, and t — s, re- 
spectively. Usually x(4) is obtained quicker from 11.13 than from 11.12, 
since the former equation involves one less matrix multiplication. How- 
ever, the integrals arising in 11.13 are generally more difficult to evalu- 
ate than those in 11.12. 


Solution with No Initial Conditions 
If no initial conditions are prescribed, the solution of X(t) = Ax(f) + f(A) is 
x(t)=e"k+ eM Te“ t(n)dt (11.15) 


or, when f(1) = 0, 
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x(t)=e'k (11.16) 
where k is an arbitrary constant vector. All constants of integration can 


be disregarded when computing the integral in Equation 11.15, since they 
are already included in k. 


Solved Problems 
Solved Problem 11.1 Put the initial-value problem 
X+2x-8x=0; x(l)=2, x(l)=3 


into the form of System 11.3. 

Following Step 1, we write ¥ = —2% + 8x; hence a,(t) =—2, a,(t) = 8, and 
f() = 0. Then, defining x,(¢) = x and x,(1) = x (the differential equation is 
second-order, so we need two new variables), we obtain x, = x,. Follow- 
ing Step 3, we find 


1 


Thus, 
Xx, = Ox, + Lx, 


Xq = 8x, — 2X 


These equations are equivalent to the matrix equation x(f) = A(x(A) + f() 


if we define 
_ x, (1) Aion: 0 1 ie 0 
x(t) = asf) (= G28 (= 0 


The differential equation is then equivalent to the matrix equation x(f) = 
A(f)x(t) + f(O), or simply x(f) = A(Dx(4), since f(A) = 0. The initial condi- 
tions can be given by x(f,) = ¢, if we define ¢, = 1 and ¢= 3 
Solved Problem 11.2 Solve x + 2x — 8x =0; x(1) = 2, x(1) =3. 
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From Problem 11.1, this initial-value problem is equivalent to Equation 
11.11 with 


SPOT eNO TA icra ok td 
x=] | AMFlg | fO=0 e=]3] = 


The solution to this system is given by Equation 11.14. For this A, e4’ is 
given in Problem 10.4; hence, 


2-1 -A(t-1 2(t-1 —A(t-1 
eAlito) — 2At-I) me 4e2@D 4 ae AED 92D _ gAe-D 
6] 8e20-D _ge-4E-D —9920-D 4 ge AU-D 
Therefore, 
x(SerrMe 


1 pam 42e 4D gD _ 4-1) ia 


“6 3 


1 2(4e24-D pe De 4 3(e°ED = etd) | 


Be2t-D _ got) 992-1) 4 got“) 


er 28624) — BoD) 4. 3(2620-) 4. de HY) 


11 574,- 1 Lage 
+f g2O-D A(t-1) 


—e 
=o ¢ 
FF 20-1) _ 2 4G-D 
6 6 


and the solution to the original initial-value problem is 


UW oq, Lacy 
x(t) = x,(t)=—e +e 
(t) = x(t) 5 6 
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Second-Order Linear Equations 
with Variable Coefficients 


A second-order linear differential equation 
by (x)y" +B (x)y’ + bo (~)y = g(x) (12.1) 


has variable coefficients when b,(x), b,(x), and bo (x) are not all constants 
or constant multiples of one another. If b,(x) is not zero in a given inter- 
val, then we can divide by it and rewrite Equation 12.1 as 


y” + P(x)y’+ Ox)y = ox) (12.2) 


where P(x) = b,(x)/b,(x), Q(x) = by(x)/b, (x), and (x) = 9(x)/b,(x). In this 
chapter, we describe procedures for solving many equations in the form 
of 12.1 and 12.2. These procedures can be generalized in a straightfor- 
ward manner to solve higher-order linear differential equations with vari- 
able coefficients. 


Analytic Functions and Ordinary Points 


A function f(x) is analytic at x, if its Taylor series about xo, 


y fF (Xp (x- 2X9)" 


n} 


n=0 


converges to f(x) in some neighborhood of x,. 


You Need to Know! 


Polynomials, sin x, cos x, and e* are analytic every- 
where. 


Sums, differences, and products of polynomials, sin x, cos x, and e* are 
also analytic everywhere. Quotients of any two of these functions are an- 
alytic at all points where the denominator is not zero. 
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The point x, is an ordinary point of the differential equation 12.2 if 
both P(x) and Q(x) are analytic at x,. If either of these functions is not an- 
alytic at x,, then x, is a singular point of 12.2. 


Solutions Around the Origin of Homogeneous 
Equations 


Equation 12.1 is homogeneous when g(x) = 0, in which case Equation 
12.2 specializes to 


y” + P(x)y’ + O(x)y =0 (12.3) 


Theorem 12.1. If x = 0 is an ordinary point of Equation 12.3, then the 
general solution in an interval containing this point has the form 


y= Y)ayx" = ayy (0) + ayy2(2) (12.4) 
n=0 
where a, and a, are arbitrary constants and y,(x) and y,(x) are linearly in- 
dependent functions analytic at x = 0. 


To evaluate the coefficients a, in the solution furnished by Theorem 12.1, 
use the following five-step procedure known as the power series method. 


Step 1. Substitute into the left side of the homogeneous differential 
equation the power series 


y= yaa = dy tajxtayx’ +a3,x° +4,x4 ++ 
o (12.5) 


n n+l n+2 
+a,Xx + Ay X + AyyoX +::- 


together with the power series for 


Y= at 2a,xt+3a;,x? +4ayx? +e 
n-l n a (12.6) 
tna,x" +(n4+layyjX" +(N4+2)dyyyx $e 


and 
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y= 2a, +6a,x+ 12a,x* eee 
+n(n—la, x"? +(ntI(n)a,,)x" (12.7) 
+(n+2)(nt Vayyyx" +++ 


Step 2. Collect powers of x and set the coefficients of each power of 
x equal to zero. 

Step 3. The equation obtained by setting the coefficient of x” to zero 
in Step 2 will contain a. terms for a finite number of j values. Solve 
this equation for the a, term having the largest subscript. The re- 
sulting equation is known as the recurrence formula for the given 
differential equation. 

Step 4. Use the recurrence formula to sequentially determine a j= 
2,3,4,...) in terms of a, and a,. 

Step 5. Substitute the coefficients determined in Step 4 into Equa- 
tion 12.5 and rewrite the solution in the form of Equation 12.4. 


The power series method is only applica- 
ble when x = 0 is an ordinary point. Although 
the differential equation must be in the form 
of Equation 12.2 to determine whether x = 0 is 
an ordinary point, once this condition is veri- 
fied, the power series method can be used on 
either form 12.1 or 12.2. If P(x) or Q(x) in 12.2 
are quotients of polynomials, it is often sim- 
pler first to multiply through by the lowest 
common denominator, thereby clearing frac- 
tions, and then to apply the power series method to the resulting equation 
in the form of Equation 12.1. 


Solutions Around the Origin of 
Nonhomogeneous Equations 


If @() in Equation 12.2 is analytic at x = 0, it has a Taylor series expan- 
sion around that point and the power series method given above can be 
modified to solve either Equation 12.1 or 12.2. In Step 1, Equations 12.5 
through 12.7 are substituted into the left side of the nonhomogeneous 
equation; the right side is written as a Taylor series around the origin. 
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Steps 2 and 3 change so that the coefficients of each power of x on the 
left side of the equation resulting from Step 1 are set equal to their coun- 
terparts on the right side of that equation. The form of the solution in Step 
5 becomes 


Y= dy, (X) + ayo (x) + y3 (xX) 


which has the form specified in Theorem 4.4. The first two terms com- 
prise the general solution to the associated homogeneous differential 
equation while the last function is a particular solution to the nonhomo- 
geneous equation. 


Initial-Value Problems 


Important! 


Solutions to initial-value problems are obtained by 
first solving the given differential equation and then 
applying the specified initial conditions. 


Solutions Around Other Points 


When solutions are required around the ordinary point x, # 0, it general- 
ly simplifies the algebra if x, is translated to the origin by the change of 
variables t= x — x9. The solution of the new differential equation that re- 
sults can be obtained by the power series method about f = 0. Then the 
solution of the original equation is easily gotten by back-substitution. 


Regular Singular Points 


The point x, is a regular singular point of the second-order homogeneous 
linear differential equation 


y” + P(x)y’ + O(x)y =0 (12.8) 
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if xj is not an ordinary point but both (x — x))P(x) and (x — Xo) O(x) are 
analytic at x,. We only consider regular singular points at x,= 0; if this is 
not the case, then the change of variables t = x — x, will translate x, to the 
origin. 


Method of Frobenius 


Theorem 12.2. If x= 0 is a regular singular point of 12.8, then the equa- 
tion has at least one solution of the form 


_ A n 
yH=x De: 
n=0 


where A and a,,(n = 0,1,2,...) are constants. This solution is valid in an in- 
terval 0 <x < R for some real number R. 


To evaluate the coefficients a, and A in Theorem 12.2, one proceeds as in 
the power series method described above. The infinite series 


with its derivatives 


y= Rig +(A+ Da,x* +(A+ Das tees 
+(A+n)a,x"" | (12.10) 


+... 


+(Atn-la,_ x4"? 


+(Atntl)a, 4.x" 


and 


y” =A(A—Dagx*? +(A4I(A)ayx* | +(A42)(A 4 Dayx" 


tet(Atn-I(Atn—2)a,_ x"? (12.11) 


+HAtn\(Atn—la, xt"? + (At nt (Atn)ay xr | t+ 
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are substituted into Equation 12.8. Terms with like powers of x are col- 
lected together and set equal to zero. When this is done for x” the result- 
ing equation is a recurrence formula. A quadratic equation in A, called the 
indicial equation, arises when the coefficient of x° is set to zero and dy is 
left arbitrary. 

The two roots of the indicial equation can be real or complex. If com- 
plex they will occur in a conjugate pair and the complex solutions that 
they produce can be combined (by using Euler’s relations and the identi- 
ty x7 = x4e*™* ) to form real solutions. In this book we shall, for 
simplicity, suppose that both roots of the indicial equation are real. Then, 
if A is taken as the /arger indicial root, A= 4, 2 ,, the method of Frobe- 
nius always yields a solution 


yaa" Ya, (A,)x" (12.12) 
n=0 
to Equation 12.8. [We have written a_(A,) to indicate the coefficients pro- 
duced by the method when 2 = ,.] 
If P(x) and Q(x) are quotients of polynomials, it is usually easier first 
to multiply 12.8 by their lowest common denominator and then to apply 
the method of Frobenius to the resulting equation. 


General Solution 


The method of Frobenius always yields one solution to Equation 12.8 of 
the form 12.12. The general solution (see Theorem 4.2) has the form y = 
C,¥,(&) +.c,y5(x) where c, and c, are arbitrary constants and y,(x) is a sec- 
ond solution of 12.8 that is linearly independent from y, (x). The method 
for obtaining this second solution depends on the relationship between 
the two roots of the indicial equation. 

Case 1. If A, —A,, is not an integer, then 


(=x? 4,0)" (12.13) 
n=0 
where y,(x) is obtained in an identical manner as y,(x) by the method 
of Frobenius, using /, in place of A,. 
Case 2. If, 2, =A,, then 


yo (x) = y, (x) Inx + x” x b, (Ay)x” (12.14) 
n=0 
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To generate this solution, keep the recurrence formula in terms of 
and use it to find the coefficients a, (n = 1) in terms of both / and a,, 
where the coefficient a, remains arbitrary. Substitute these a, into 
Equation 12.9 to obtain a function y(A, x) which depends on the vari- 
ables A and x. Then 


dy(A, x) 
x)= 12.15 
(x) OA ie ( ) 
Case 3. If, —A, =N, a positive integer, then 
yp (x) = d_yy,(x)Inx +x” ¥ d,(Ay)x” (12.16) 


n=0 


To generate this solution, first try the method of Frobenius with /,. 
If it yields a second solution, then this solution is y,(x), having the 
form of 12.16 with d_, =0. Otherwise, proceed as in Case 2 to gen- 
erate y(A, x), whence 


(x)= S1- AA, (12.17) 


Solved Problems 


Solved Problem 12.1 Determine whether x = 0 is an ordinary point of 
the differential equation 


y”—xy’+2y=0 


Here P(x) = —x and Q(x) = 2 are both polynomials; hence they are 
analytic everywhere. Therefore, every value of x, in particular x = 0, is 
an ordinary point. 


Solved Problem 12.2 Find a recurrence formula for the power series so- 
lution around x = 0 for the differential equation given in Problem 12.1. 


It follows from Problem 12.1 that x = 0 is an ordinary point of the given 
equation, so Theorem 12.1 holds. Substituting Equations 12.5 through 
12.7 into the left side of the differential equation, we find 
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[2a, + 6a,x + Dax +e $n(n— has +(nt+ 1)(n)a, 4.x" tee] 
—x[a, + 2a)x+3a5x7 +++--+na, x") +(n+Na, yx" +] 
+2[dy + ayx+ yx? +03x? +---44,x" +a, x" ++] =0 


Combining terms that contain like powers of x, we have 


(2d, + 2ay)+ x(6a, +a,)+ x” (12a,)+x° (20a, — a3) 
te + x"[(nt2)(n+ Da, 45 —na, +2a,]+-- 


=0+0x+0x7 40x73 +--+ 0x" 4-5 


The last equation holds if and only if each coefficient in the left-hand side 
is zero. Thus, 


2a, + 2a =0, 6a3 +a =0, 12a, =0, 20as5 — a3 =0, wine 
In general, 


(n+2)(n+)a,,.—-(n—-2)a, =0, or, 
(n—2) 
Qn+2 = avn 
(n+2)(n+1) 


which is the recurrence formula for this problem. 


Solved Problem 12.3 Find the general solution near x = 0 of y" — xy’ + 
2y=0. 


Successively evaluating the recurrence formula obtained in Problem 12.2 
for n = 0,1,2,..., we calculate 


ay =—ay 
1 
a, Pet 
a4=0 (12.18) 
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2 
Ege 

M5 = 39 4 = 75) 

a= 3 ac = al u a, |= u a 

™ 42> 144 120 |) 1680 ! 
4 1 

= =—(0)=0 
Br se ig 


Note that since a, = 0, it follows from the recurrence formula that all the 
even coefficients beyond a, are also zero. Substituting 12.18 into Equa- 
tion 12.5 we have 


2 3 4 5 6 7 
=d),+a,x-d)x° —-~a,x° +0x" —-—~a,x° +0x° - 
pe 6 | 120 | 


=ay(l- 3°) +(x ax wor bf --) 


x 
120 1680 


If we define 


Bath 5 1 Pe 
120 1680 


y(ay=l-x? and y,(x)=x a 


then the general solution can be rewritten as y = ayy, (x) + a,y5(x). 


Solved Problem 12.4 Determine whether x=0 is a regular singular point 
of the differential equation 


2.” 


8x°y" +10xy’+(x-Dy=0 


Dividing by 8x”, we have 


2) 1 1 
P(x)=— and x)=—-—> 
() 4x Qe) 8x 8x? 
Neither of these functions is defined at x = 0, so this point is a singular 
point. Furthermore, both 


_> eee 
aECO= and HOGS 1) 
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are analytic everywhere: the first is a constant and the second a polyno- 
mial. Hence, both are analytic at x = 0, and this point is a regular singu- 
lar point. 


Solved Problem 12.5 Use the method of Frobenius to find one solution 
near x =0 of x°y”—xy’+y=0. 


Here P(x) =—1/x and Q(x) = 1/x?, so x = 0 is a regular singular point 
and the method of Frobenius is applicable. Substituting Equations 12.9 
through 12.11 into the left side of the differential equation, as given and 
combining coefficients of like powers of x, we obtain 


x*(A-1) ay + x4" [V7 a,]+-- 
+x"*"[(A+n) —2(A+n)+l1la, +:-=0 
Thus, 
(A-1)' a, =0 (12.19) 
and, in general, 
[(At+n)* —2+n)+1]a, =0 (12.20) 
From Equation 12.19, the indicial equation is (A — 1)? = 0, which has roots 
A, =A, = 1. Substituting 2 = 1 into Equation 12.20 we obtain na, = 0, 


which implies that a, = 0, n 2 1. Thus, y, (x) = apx. 


Solved Problem 12.6 Find the general solution near x = 0 of 3x*y”—- 
xy’+y=0. 

Here P(x) =—1/(3x) and Q(x) = 1/(3x?); hence, x = 0, is a regular sin- 
gular point and the method of Frobenius is applicable. Substituting Equa- 


tions 12.9 through 12.11 into the differential equation and simplifying, 
we have 


x*[3A* — 424 lay +.x4"'[327 + 2A]a, ++ 
+x"*"[3(A +n) —4(A+n)+1]a, +---=0 


Dividing by x? and equating all coefficients to zero, we find 
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(3A? —4A+1)ay =0 (12.21) 


and 


[3A+n)—4A+n)+lJa,=0 n21 (12.22) 


From 12.21, we conclude that the indicial equation is 3A? — 42 + 1=0, 
which has roots 2, = 1 and 2, = 1/3. Since A, — 2, = 2/3, the solution is 
given by Equations 12.12 and 12.13. Note that for either value of A, Equa- 
tion 12.22 is satisfied by simply choosing a, = 0,n 2 1. Thus, 


y= x 3 A,X" = AyX 
n=0 


173 173 
y7(x) =x > A,X" = Agx 
n=0 


and the general solution is 
y=C1y1 (4) + Oy) = hx t yx" 


( 


where ky =C\d and k, =C dy 


Solved Problem 12.7 Use the method of Frobenius to find one solution 


near x=0 of x*y”+(x* —2x)y’+2y=0. 


Here P(x)= ie and Q(x)= ae so, x = 0 is a regular singular 
x x 


point and the method of Frobenius is applicable. Substituting Equations 
12.9 through 12.11 into the left side of the differential equation, as giv- 
en, and combining coefficients of like powers of x, we obtain 


x*((A? —3A + 2)ag]+x*"[(A? — Aja, + ag] +- 
+x" ([(Atn) —3(Atn)+2]a,+(Atn—-Da,_}+-=0 


Dividing by x*, factoring the coefficient of a,, and equating the coeffi- 
cient of each power of x to zero, we obtain 


(A? —3A+2)ay =0 (12.23) 
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and, in general, [((A+n)—2][(A+n)—-l]a, +(At+n—Da,_, =9, or, 


1 
a, =->——4,-, (n21 12.24 
pee ay ee ( ) ( ) 
From 12.23, the indicial equation is A? — 3A + 2 = 0, which has roots A, 
=2and/, = 1. Since A, —2,= 1, a positive integer, the solution is given 
by Equations 12.12 and 12.16. Substituting 2 = 2 into 12.24, we have a, 


=—(1/n)a,_,, from which we obtain 
a, =—ao 
a —— = 
OE Oe 
—_ 1 _ ol By = 1 
nce Bees. gDree cape 
: Cf 
and, in general, a, = dy. Thus, 


k 


— (-1)" 2 
yx) = ax? dS ) x” =ayx’e* 
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Series 

Y Solved Problems 


Gamma Function 
The gamma function, T(p), is defined for any positive real number p by 
Pp) = fx ea (13.1) 
Consequently, (1) = 1 and for a positive real number p, 
(p+ D= p(y) (13.2) 
Furthermore, when p = n, a positive integer, 


T(n+1)=n! (13.3) 


98 
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Thus, the gamma function (which is defined on all positive real numbers) 
is an extension of the factorial function (which is defined only on the non- 
negative integers). 

Equation 13.2 may be rewritten as 


1 
HAP) 5 SB ae (13.4) 


which defines the gamma function iteratively for all nonintegral negative 
values of p. [(0) remains undefined, because 


ony S te 
p-0* p>0* 
and 
indie See 
p00” po0- p 


It then follows from Equation 13.4 that I’(p) is undefined for negative in- 
teger values of p. 

Table 13.1 lists values of the gamma function in the interval 1 < p< 
2. These tabular values are used with Equations 13.2 and 13.4 to gener- 
ate values of I(p) in other intervals. 


Bessel Functions 


Let p represent any real number. The Bessel function of the first kind of 
order p, J), is 


5 (—1)* x2k+P 
J,(x) = (13.5) 
z A 2? PIT (p+k +) 


The function JX) is a solution near the regular singular point x = 0 of 
Bessel’s differential equation of order p: 


xy” + xy’ +(x? — p?)y=0 (13.6) 


In fact, JX) is that solution of Equation 13.6 guaranteed by Theorem 
12.2. 
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— 


x TQ) x T(x) x T(x) 
1.00 | 1.0000 0000 | 1.34 | 0.8922 1551 | 1.67 | 0.9032 9650 
1.01 | 0.9943 2585 | 1.35 | 0.89115144 | 1.68 | 0.9050 0103 
1.02 | 0.9888 4420 | 1.36 | 0.8901 8453 | 1.69 | 0.9067 8182 
1.03 | 0.9835 4995 | 1.37 | 0.8893 1351 | 1.70 | 0.9086 3873 
1.04 | 0.9784 3820 | 1.38 | 0.8885 3715 | 1.71 | 0.9105 7168 
1.05 | 0.9735 0427 | 1.39 | 0.88785429 | 1.72 | 0.9125 8058 
1.06 | 0.9687 4365 | 1.40 | 0.8872 6382 | 1.73 | 0.9146 6537 
1.07 | 0.9641 5204 | 1.41 | 0.8867 6466 | 1.74 | 0.9168 2603 
1.08 | 0.9597 2531 | 1.42 | 0.88635579 | 1.75 | 0.9190 6253 
1.09 | 0.9554 5949 | 1.43 | 0.8860 3624 | 1.76 | 0.9213 7488 
1.10 | 0.9513 5077 | 1.44 | 0.8858 0506 | 1.77 | 0.9237 6313 
1.11 | 0.9473 9550 | 1.45 | 0.8856 6138 | 1.78 | 0.9262 2731 
1.12 | 0.9435 9019 | 1.46 | 0.8856 0434 | 1.79 | 0.9287 6749 
1.13 | 0.9399 3145 | 1.47 | 0.8856 3312 | 1.80 | 0.9313 8377 
1.14 | 0.9364 1607 | 1.48 | 0.8857 4696 | 1.81 | 0.9340 7626 
1.15 | 0.9330 4093 | 1.49 | 0.8859 4513 | 1.82 | 0.9368 4508 
1.16 } 0.9298 0307 ; 1.50 | 0.8862 2693 | 1.83 | 0.9396 9040 
1.17 | 0.9266 9961 | 1.51 | 0.8865 9169 | 1.84 | 0.9426 1236 
1.18 | 0.9237 2781 ; 1.52 | 0.88703878 | 1.85 | 0.9456 1118 
1.19 | 0.9208 8504 | 1.53 | 0.8875 6763 | 1.86 | 0.9486 8704 
1.20 | 0.9181 6874 | 1.54 | 0.8881 7766 | 1.87 | 0.9518 4019 
1.21 | 0.9155 7649 | 1.55 | 0.8888 6835 | 1.88 | 0.9550 7085 
1.22 | 0.9131 0595 | 1.56 | 0.8896 3920 | 1.89 | 0.9583 7931 
1.23 | 0.9107 5486 | 1.57 | 0.8904 8975 | 1.90 | 0.9617 6583 
1.24 | 0.9085 2106 | 1.58 | 0.89141955 | 1.91 | 0.9652 3073 
1.25 | 0.9064 0248 | 1.59 | 0.8924 2821 | 1.92 | 0.9687 7431 
1.26 | 0.9043 9712 | 1.60 | 0.8935 1535 {| 1.93 | 0.9723 9692 
1.27 | 0.9025 0306 | 1.61 | 0.8946 8061 | 1.94 | 0.9760 9891 
1.28 | 0.9007 1848 | 1.62 | 0.8959 2367 | 1.95 | 0.9798 8065 
1.29 | 0.8990 4159 | 1.63 | 0.8972 4423 | 1.96 | 0.9837 4254 
1.30 | 0.8974 7070 | 1.64 | 0.8986 4203 | 1.97 | 0.9876 8498 
1.31 | 0.8960 0418 | 1.65 | 0.9001 1682 | 1.98 | 0.9917 0841 
1.32 | 0.8946 4046 | 1.66 | 0.9016 6837 | 1.99 | 0.9958 1326 
1.33 | 0.8933 7805 


Table 13.1 
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Algebraic Operations on Infinite Series 


Changing the dummy index. The dummy index in an infinite series can 
be changed at will without altering the series. For example, 


— — 1 — 1 1 1 1 #1 ~=41 
py = =—+—+—4+—4 4... 
Pan cei Aviat)! AA(ptD! Wo 2 Bt 4 St 


1 
Change of variables. Consider the infinite series eer 4 k+ pr - If we make 


the change of variables j=k+1, ork=j—1, then “ 


co 


— 1 
X Ee 1)! Xi ‘jl 
Note that a change of variables generally changes the limits on the sum- 
mation. For instance, if j =k + 1, it follows that j = 1 when k = 0, j =<, 
when k = 9, and, as k runs from 0 to ©, j runs from | to ©, 

The two operations given above are often used in concert. For ex- 
ample, 


co 


py rm ap oy 7 ey Ge = 


jr2 k=2 


Here, the second series results from the change of variables j = k + 2 in 
the first series, while the third series is the result of simply changing the 
dummy index in the second series from j to k. Note that all three series 
equal 


Solved Problems 
Solved Problem 13.1 Prove that [(p + 1) =pI(p), p > 0. 


Using Equation 13.1 and integration by parts, we have 
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T(p+l= jx (pt)-1, “ae tim fx e “dx 


reo 


= lim| —x?e* e “dx 


rod 


= lim(-r?e" +0)+ pj x? le“ dx = pI(p) 
roo 
0 
The result lim r?e" 
roo 
and using L’ H6pital’s rule. 


= 0 is easily obtained by first writing r?e~" as r?/e” 


Solved Problem 13.2 Use the method of Frobenius to find one solution 
of Bessel’s equation of order p: 


xy" + xy’ +(x? — p”)y =0 


Substituting Equations 12.9 through 12.11 into Bessel’s equation and 
simplifying, we find that 


x4 (A? — p”)ay +x" [A +1) — p? Ja, + 


$x? ([(A +2)" =p? Jay tay} +e 
+x" ((A+n) — pla, +4,-5} + =0 
Thus, 
(A? =p? )ay =0 [(A+1)? —p? Ja, =0 (13.7) 


and, in general, [(A+ ny - pla, +d,-7 =0, or 


1 
a, =-———— yn (N22 13.8 
GamPapeine (22) (13.8) 
The indicial equation is A? — p?=0, which has roots 4, = p and 2, =—p 
(p nonnegative). Substituting 2 = p into 13.7 and 13.8 and simplifying, 
we find that a, = 0 and 


a a es (n2 2) 
n(2p+n) 
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Hence, 0=a, =a; =a; =a, =--: and 
ol 
Oo DA p+h © 
ay =-a : = 4 : A 
2°2(pt2) 2°2!(p+2)(p+1) 
ne 1 -1 


P3(p+3) ' 23\(p+3(p+D(p+) ° 


and, in general, 


—1)* 
ee) ee 
2° kK\(pt+k)(ptk—-1)---(p+2)(pt+D 
y,(x) =x? by a,x" = ake + Soc 
n=0 k=1 
Thus, 
co yk 2k 
=x? i+) xT (oD 
m2. KM ptk)\(pt+k—l)---(p+2)(pt+) 
1 
It is customary to choose the arbitrary constant a, as ay = ——————. 
3 z Va aa A Dee ah 


Then bringing a,x? inside the brackets and summation, combining, and 
finally using Problem 13.1, we obtain 
= (—1)k x24? 


1 
xP + 
2°T(p+1) LF rr sks 


y= 


oo ei 2k+P : 
-> aia = J,(x) 
te T(pt+k+)) 
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Direction Fields 


Graphical methods produce plots of solutions to first-order differential 
equations of the form 


y = f(x,y) (14.1) 
where the derivative appears only on the left side of the equation. 


Example 14.1 (a) For the problem y’ = —y + x + 2, we have f(x, y) =—y 
+x +2. (b) For the problem y' = y + 1, we have f(x, y) = y* + 1. (c) For 
the problem y’ = 3, we have f(x, y) =3. Observe that in a particular prob- 
lem, f(x, y) may be independent of x, of y, or of x and y. 


Equation 14.1 defines the slope of the solution curve y(x) at any point 
(x, y) in the plane. A line element is a short line segment that begins at the 
point (x, y) and has a slope specified by 14.1; it represents an approxi- 
mation to the solution curve through that point. A collection of line ele- 
ments is a direction field. Note that direction fields are tedious to draw by 
hand, but they can be easily created by many computer algebra systems. 
The graphs of solutions to 14.1 are generated from direction fields by 
drawing curves that pass through the points at which line elements are 
drawn and also are tangent to those line elements. 


wx Note! 


A direction field can be used to graphically deter- 
mine the behavior of the solution. 


If the left side of Equation 14.1 is set equal to a constant, the graph of the 
resulting equation is called an isocline. Different constants define differ- 
ent isoclines, and each has the property that all elements emanating from 
points on that isocline have the same slope, a slope equal to the constant 
that generated the isocline. When they are simple to draw, isoclines yield 
many line elements at once which are useful for constructing direction 
fields. 
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Euler’s Method 


If an initial condition of the form 


W(Xo) = Yo (14.2) 


is also specified, then the only solution curve of Equation 14.1 of inter- 
est is the one that passes through the initial point (x), yo). 

To obtain a graphical approximation to 
the solution curve of Equations 14.1 and 
14.2, begin by constructing a line element at 
the initial point (x, y,) and then continuing 
it for a short distance. Denote the terminal 
point of this line element as (x,, y,). Then 
construct a second line element at (x,, y,) and 
continue it a short distance. Denote the ter- 
minal point of this second line element as (x,, 
y,). Follow with a third line element con- 
structed at (x,, y,) and continue it a short distance. The process proceeds 
iteratively and concludes when enough of the solution curve has been 
drawn to meet the needs of those concerned with the problem. 

If the difference between successive x values are equal, that is, if for 
a specified constant h, h= x, — x) =x, —x, =x, —x, =-+:, then the graph- 
ical method given above for a first-order initial-value problem is known 
as Euler’s method. It satisfies the formula 


Yas =Yn thf (Xn Yn) (14.3) 


for n = 1,2,3,.... This formula is often written as 


Yntl =Vn thy, (14.4) 


where 


Yn = SX Yn) (14.5) 


as required by Equation 14.1. 


CHAPTER 14: Numerical Methods 107 
Stability 


The constant h in Equations 14.3 and 14.4 is called the step-size, and its 
value is arbitrary. In general, the smaller the step-size, the more accurate 
the approximate solution becomes at the price of more work to obtain that 
solution. Thus, the final choice of h may be a compromise between ac- 
curacy and effort. If is chosen too large, then the approximate solution 
may not resemble the real solution at all, a condition known as numeri- 
cal instability. To avoid numerical instability, Euler’s method is repeat- 
ed, each time with a step-size one-half its previous value, until two suc- 
cessive approximations are close enough to each other to satisfy the needs 
of the solver. 


General Remarks Regarding 
Numerical Methods 


A numerical method for solving an initial-value problem is a procedure 
that produces approximate solutions at particular points using only the 
operations of addition, subtraction, multiplication, division, and func- 
tional evaluations. In this chapter, we consider only first-order initial-val- 
ue problems of the form 


y =flxy); yx) = Yo (14.6) 


Generalizations to higher-order problems are given later in this chapter. 
Each numerical method will produce approximate solutions at the points 
Xp»X1,%,-.-, Where the difference between any two successive x-values is 
a constant step-size h; that is x,,, —x, =A (n=0,]1,2,...). Remarks made 
previously in this chapter on the step-size remain valid for all the nu- 
merical methods presented. 

The approximate solution at x, will be designated by y(x,), or sim- 
ply y,,. The true solution at x, will be denoted by either Y(x,,) or Y_. Note 
that once y_ is known, Equation 14.6 can be used to obtain y’ as 


Yn =f (%nsYn) (14.7) 


The simplest numerical method is Euler’s method, described earlier 
in this chapter. 
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A predictor-corrector method is a set of two equations for y,_,. The 
first equation, called the predictor, is used to predict (obtain a first ap- 
proximation to) y,,,,; the second equation, called the corrector, is then 
used to obtain a corrected value (second approximation) to y,,,. In gen- 


eral, the corrector depends on the predicted value. 


Modified Euler’s Method 


This is a simple predictor-corrector method that uses Euler’s method as 
the predictor and then uses the average value of y’ at both the left and right 
end points of the interval [x,, x, ,,] (2 =0,1,2,...) as the slope of the line 


wor 
element approximation to the solution over that interval. The resulting 
equations are: 


predictor: y,,, =y, thy, 


h , tf 
corrector: yn.) =), + 5 Ont ty 


For notational convenience, we designate the predicted value of y,,,, by 
Py,,,;- It then follows from Equation 14.7 that 


PY n+ = f (Xn41> PYnsi) (14.8) 
The modified Euler’s method becomes 


predictor: py,4)=y, thy, 


corrector: y,4)=y,+ + (Ya he y,) (14.9) 
Runge-Kutta Method 
Yat =In te (ky + 2k, + 2k; +k,) (14.10) 


where 


ky =hf (Xn.Yn) 
1 1 
k, =hf(x, thy, res 


2 
ky =hf(x, thy, +k) 
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This is not a predictor-corrector method. 


Adams-Bashforth-Moulton Method 


predictor: py,.,=y,+ a (55y, —59y"_,+37y;_5 —9y7_3) 


i‘ (14.11) 
corrector: Yj.) =Y,+ 54 OPYntt +19y, — 5Yn—1 + Yn-2) 
Milne’s Method 
. 4h , vy tl 
predictor: py,4) =Y,-3+ 3 (2), 7 Yat 2) 2) 
(14.12) 


h , , , 
corrector: y,4)=Y,-)+ 3 (PYnsi + 4Vn + Yn-1) 
Starting Values 


The Adams-Bashforth-Moulton method and Milne’s method require in- 
formation at yp, y,, Y>, and y, to start. The first of these values is given by 
the initial condition in Equation 14.6. The other three starting values are 
gotten by the Runge-Kutta method. 


Order of a Numerical Method 


A numerical method is of order n, where n is a positive 

integer, if the method is exact for polynomials of de- © 
gree n or less. In other words, if the true solution of an ZA 
initial-value problem is a polynomial of degree n or /< Sg 
less, then the approximate solution and the true solu- 

tion will be identical for a method of order n. 

In general, the higher the order, the more accurate the method. 
Euler’s method, Equation 14.4, is of order one, the modified Euler’s 
method, Equation 14.9, is of order two, while the other three, Equations 
14.10 through 14.12, are fourth-order methods. 


\ 
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Numerical Methods for Systems 


First-Order Systems 


Numerical methods for solving first-order initial-value problems, includ- 
ing all of those described previously in this chapter, are easily extended 
to a system of first-order initial-value problems. These methods are also 
applicable to most higher-order initial-value problems, in particular 
those that can be transformed to a system of first-order differential equa- 
tions by the reduction process described in Chapter Eleven. 

Standard form for a system of two equations is 


y= f(x, y,2) 
Zz’ = (x,y,z); (14.13) 
WXx9)= Yo» 2%) = Zo 
We note that, with y’ = f(x, y,z)=z, with system 14.13 represents the 
second-order initial-value problem 


y”=e(% yy); W%XO)=Yo, (Xo) = Zo 


Standard form for a system of three equations is 


y= f(x, y,zw) 
Zz’ = g(x, y,Z,W) 
w= r(x, y,Z,W); (14.14) 


V(X )=Yo» WX) =Zy, W(X) = Wo 


If, in such a system, f(x, y, Zz, w) = z, and g(x, y, z, w) = w, then system 
14.14 represents the third-order initial-value problem 


y” = r(x, y,Z,W)} y(X9) = Yo. y’ (Xo) = Z05 y’ (Xp) = Wo 


The formulas that follow are for systems of two equations in stan- 
dard form (14.13). Generalizations to systems of three equations in stan- 
dard form (14.14) or systems with four or more equations are straight- 
forward. 


CHAPTER 14: Numerical Methods 111 


Euler’s Method 


=y, thy; 
Yn+t = Yn Ms (14.15) 
Zn+l = Sn +hz, 


Runge-Kutta Method 


Yast =n +o tk + 2k, +2k3 + ky) 


1 (14.16) 
Zntl = Zn ee +21, +21, +1,) 


where 


Ky = Pf (Xn YnZn) 

Ly =Ng(Xn»Vno2n) 

ky = hf (X_ + Zn t+ Then + Zh) 
by =hg(Xn t+ Zhen t+ 7 hln + Zh) 
ky =hf (i, + hn t+ then t Fh) 
I, =hg(x, +51, y, + oho kn toh) 
ky =hf(x, thy, +k3,Z, +) 

I, = hg(x, thy, +k,z, +45) 


Adams-Bashforth-Moulton Method 


predictors: Yn.) =, + = (S59, = 99), +37 9 99-9) 

P2na1 =2n + s(55c; —592)_,+37z,_, -927_3) 
correctors: ny) =Yq + SPY S199 Sy yg): ED 
Zntl = %_ + = Opi 19s) 95, 4S 8) 


Corresponding derivatives are calculated from system 14.13. In particu- 
lar, 
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Yntl a ne een ere pean p) (14.18) 
Zn = 8Cns VntiZn41) 


The derivatives associated with the predicted values are obtained simi- 
larly, by replacing y and z in Equation 14.18 with py and pz, respective- 
ly. As in the previous section, four sets of starting values are required for 
the Adams-Bashforth-Moulton method. The first set comes directly from 
the initial conditions; the other three sets are obtained from the Runge- 
Kutta method. 


Solved Problems 


Solved Problem 14.1 Draw two solution curves to the differential equa- 
tion y =2y-x. 


A direction field for this equation is given by Figure 14-1. Two solution 
curves are also shown, one that passes through the point (0,0) and a sec- 
ond that passes through the point (0,2). Observe that each solution curve 
follows the flow of the line elements in the direction field. 


| | / / / / 
/ / / / / 4 
/ / / / 4 2 
/ / 7 2 _— 
ae Be pa tee 2s 
/ / Seat Se Ke % 
os 7% Se OM CARS AS 
ue Me YN ONY 
= ae ar 4 
\ \ \ \ 
\ \ \ \ 
\ \ \ \ 
\ \ \ \ 


Figure 14-1 
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Solved Problem 14.2 Find y(1) for y’=y—.x; y(O) = 2, using Euler’s 
method with h=+4 


For this problem, x, = 0, yy = 2, and f(x, y) = y — x; so Equation 14.5 be- 
comes y;, =y, —x,- Because h= 7, 


Using Equation 14.4 with n = 0,1,2,3 successively, we now compute the 
corresponding y-values. 


n=0: y= Yo + hy 
But y= So, Yo) = Yo X= 2-0 = 2 
Hence, y, =2+4 (2)=3 

n=1: Yo =y, thy, Se 
But y, =f@,¥,) =), -*%,= F-4=7 
Hence, y, = $+1(2)=# 

n=2: ¥3=Yy + hy, r 
Beenie) x,=2 1 - 4l 
Hence, y3=7% 244 14 1) = 237 


n=3: y= yz + hy; 
But y= F(t Ys) = Ya x,= t= =@ 


= 33 189) _ i137 
Hence, y, + aera )= a5 


Thus, 
YL) = 4 = 956 = 4.441 


Note that the true solution is Y(x) = ex+x+1,so that YC) = 4.718. 
If we plot (x,, y,,) for n = 0,1,2,3, and 4, and then connect successive 
points with straight line segments, as done in Figure 14-2, we have an ap- 
proximation to the solution curve on [0,1] for this initial-value problem. 
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Figure 14-2 


Chapter 15 
BOUNDARY- VALUE 


PROBLEMS AND 
FOURIER SERIES 


In THIS CHAPTER: 


Vv 


VC NSNLNRN 


Second-Order Boundary- Value 
Problems 

Eigenvalue Problems 
Sturm-Liouville Problems 
Eigenfunction Expansions 
Fourier Sine Series 

Fourier Cosine Series 

Solved Problems 


Second-Order Boundary-Value Problems 


Standard Form 


A boundary-value problem in standard form consists of the second-order 
linear differential equation 


y" + P(x)y’+ Oa)y = G(x) (15.1) 
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and the boundary conditions 


a,y(a)t+ Byy(ay=7, 


: (15.2) 
an y(b)+ Byy"(b) = 72 


where P(x), Q(x), and (x) are continuous in [a, b] and a,, a), B,, 8, 7;, 
and y, are all real constants. Furthermore, it is assumed that a, and f, are 
not both zero, and also that a, and £, are not both zero. 

The boundary-value problem is said to be homogeneous if both the 
differential equation and the boundary conditions are homogeneous (i.e., 
0(x) = 0 and y, = y, = 0). Otherwise the problem is nonhomogeneous. 
Thus a homogeneous boundary-value problem has the form 


y" + P(x)y’ + Q(a)y = 0; 
a y(a)+ By'(ay=0 (15.3) 
a,y(b)+ Byy'(b) =0 
A somewhat more general homogeneous boundary-value problem than 
15.3 is one where the coefficients P(x) and Q(x) also depend on an arbi- 
trary constant A. Such a problem has the form 
y" + P(x, A)y’ + O(x, A)y = 05 
a y(a)+ By’(a) =0 (15.4) 
a y(b)+ Bry'(b) = 0 
Both 15.3 and 15.4 always admit the trivial solution y(x) = 0. 
Solutions 
A boundary-value problem is solved by first obtaining the general solu- 
tion to the differential equation, using any of the appropriate methods pre- 
sented heretofore, and then applying the boundary conditions to evaluate 


the arbitrary constants. 


Theorem 15.1. Let y,(x) and y,(x) be two linearly independent solutions 
of 


y” + P(x)y’ + Q(x)y =0 
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Nontrivial solutions (i.e., solutions not identically equal to zero) to the 
homogeneous boundary-value problem 15.3 exist if and only if the de- 
terminant 


ay (a)+Byyj(a) Oyo (a)+ Bi y3(a) 


, ; (15.5) 
yy (b) + Boy (b) Oy ¥2(b)+ Boy3(d) 


equals zero. 


Theorem 15.2. The nonhomogeneous boundary-value problem defined 
by 15.1 and 15.2 has a unique solution if and only if the associated ho- 
mogeneous problem 15.3 has only the trivial solution. 


In other words, a nonhomogeneous problem has a unique solution 
when and only when the associated homogeneous problem has a unique 
solution. 


Eigenvalue Problems 


When applied to the boundary-value problem 15.4, Theorem 15.1 shows 
that nontrivial solutions may exist for certain values of 4 but not for oth- 
er values of A. 


You Need to Know J 


Those values of 4 for which nontrivial solutions do 
exist are called eigenvalues; the corresponding 
nontrivial solutions are called eigenfunctions. 


Sturm-Liouville Problems 


A second-order Sturm-Liouville problem is a homogeneous boundary- 
value problem of the form 


[p(x)y’V + q(x)y + Aw(x)y = 0; (15.6) 


118 DIFFERENTIAL EQUATIONS 


ay(a)+ Biy(a)=0 

ay(b) + Bry'(b) = 0 
where p(x), p’(x), g(x), and w(x) are continuous on [a, b], and both p(x) 
and w(x) are positive on [a, b]. 

Equation 15.6 can be written in standard form 15.4 by dividing 
through by p(x). Form 15.6, when attainable, is preferred, because Sturm- 
Liouville problems have desirable features not shared by more general 
eigenvalue problems. The second-order differential equation 


(15.7) 


ay(x)y” +.a,(x)y’ +.dg(x)y + Ar(x)y =0 (15.8) 


where a,(x) does not vanish on [a, 5], is equivalent to Equation 15.6 if 
and only if a(x) = a,(x). This condition can always be forced by multi- 
plying Equation 15.8 by a suitable factor. 


Properties of Sturm-Liouville Problems 


Property 15.1. The eigenvalues of a Sturm-Liouville problem are real 
and nonnegative. 


Property 15.2. The eigenvalues of a Sturm-Liouville problem can be 
arranged to form a strictly increasing infinite sequence; that is, 0 <A, < 


A, <A, <--:. Furthermore, 2, > o as n > o. 


Property 15.3. For each eigenvalue of a Sturm-Liouville problem, there 
exists one and only one linearly independent eigenfunction. 


[By Property 15.3 there corresponds to each eigenvalue 4,, a unique 
eigenfunction where the coefficient on the lead term is one; we denote 
this eigenfunction by e (x).] 

Property 15.4. The set of eigenfunctions {e,(x), e,(x),...} of a Sturm- 
Liouville problem satisfies the relation 


b 
woe, (x)e,,(x)dx = 0 (15.9) 


for n #m, where w(x) is given in Equation 15.6. 
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Eigenfunction Expansions 


A wide class of functions can be represented by infinite series of eigen- 
functions of a Sturm-Liouville problem. 


Definition: A function f(x) is piecewise continuous on the open interval 
a<x<bif (1) f(x) is continuous everywhere in a < x < b with the possi- 
ble exception of at most a finite number of points x,, X,,..., x,,, and (2) at 


these points of discontinuity, the right- and left-hand limits of f(x) re- 
spectively lim f(x) and lim f(x), exist (j = 1,2,..., 2). 


xX>X; X<Xj; 


(Note that a continuous function is piecewise continuous.) 


Definition: A function f(x) is piecewise continuous on the closed inter- 
val asx < bif (1) it is piecewise continuous on the open interval a<x<b, 
(2) the right-hand limit of f(x) exists at x = a and (3) the left-hand limit 
of f(x) exists at x = b. 


Definition: A function f(x) is piecewise smooth on [a, b] if both f(x) and 
f(x) are piecewise continuous on [a,b]. 


Theorem 15.3. If f(x) is piecewise smooth on [a, b] and if {e,,(x)} is the 
set of all eigenfunctions of a Sturm-Liouville problem (see Property 
15.3), then 


FOES 25H) (15.10) 
n=1 
where 
b 
J w(x) f(x)e,, (x)dx 
a (15.11) 
Jwoe; (x)dx 


The representation 15.10 is valid at all points in the open interval (a, b) 
where f(x) is continuous. The function w(x) in 15.11 is given by Equa- 
tion 15.6. 
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Because different Sturm-Liouville problems usually generate differ- 
ent sets of eigenfunctions, a given piecewise smooth function will have 


many expansions of the form 15.10. The basic features of all such ex- 
pansions are exhibited by the trigonometric series discussed below. 


Fourier Sine Series 


The eigenfunctions of the Sturm-Liouville problem y” + Ay =0; y(0) =0, 
y(L) = 0, where L is a real positive number, are e (x) = sin(nax/L) (n = 
1,2,3,...). Substituting these functions into 15.10, we obtain 


f(x)= sc, sin (15.12) 
n=1 


For this Sturm-Liouville problem, w(x) = 1, a= 0, and b= L; so that 


fnde2 (ah a ifte se = ze 
: i‘ 6 L 2 
and 15.11 becomes 
ee ae eace (15.13) 
aie if L 


The expansion 15.12 with coefficients given by 15.13 is the Fourier sine 
series for f(x) on (0, L). 


Fourier Cosine Series 


The eigenfunctions of the Sturm-Liouville problem y”’ + Ay =0; y’(0) = 
0, y(L) = 0, where L is a real positive number, are e,(x) = 1 and e (x) = 
cos(nrx/L) (n = 1,2,3,...). Here 2 = 0 is an eigenvalue with correspond- 
ing eigenfunction e,(x) = 1. Substituting these functions into 15.10, 
where because of the additional eigenfunction e,(x) the summation now 
begins at n = 0, we obtain 


F(X) = €9 + Yen C08 (15.14) 
n=1 
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For this Sturm-Liouville problem, w(x) = 1, a= 0, and b= L; so that 


b L 
woes (x)dx = Jax =L 


a 0 
and 
? ; NIX L 
Jue; (x)dx = Jcos? dx = 
‘ 5 L 2 
Thus 15.11 becomes 
ia 0 NIX 
Co=— x)dx cz x) cos — dx 15.15 
0 rt@ : races : (15.15) 
(n=1,2,...) 


The expansion 15.14 with coefficients given by 15.15 is the Fourier co- 
sine series for f(x) on (0, L). 


Solved Problems 
Solved Problem 15.1 Solve y” + 2y’ — 3y= 9x; yO) = 1, y'(1) =2. 


This is a nonhomogeneous boundary-value problem of forms 15.1 
and 15.2, where ¢(x) = 9x, y, = 1, and y, = 2. The general solution to the 
homogeneous differential equation is y = c,e~** + c,e*. If we apply ho- 
mogeneous boundary conditions, we find that c, = c, = 0; hence the so- 
lution to the homogeneous problem is y = 0. Since the associated homo- 
geneous problem has only the trivial solution, it follows from Theorem 
15.2 that the given problem has a unique solution. Solving the differential 
equation by the method of Chapter Six, we obtain 


y=ce * +ce,e* —3x-2 


Applying the boundary conditions, we find 
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Cte) -2=1 —3ce%+c,e-3=2 


whence 
3e—5 5+9e3 
3 = 3 
e+3e e+3e 


Cy 


Finally, 
y= Bex 5)e* +(5+9e™ e* 


= 3x-2 
e+3e 


Solved Problem 15.2 Find the eigenvalues and eigenfunctions of 
y"—4ay'+42*%y=0;  y(0)=0, y+ yD =0 


The coefficients of the given differential equation are constants (with 
respect to x); hence, the general solution can be found by use of the char- 
acteristic equation. We write the characteristic equation in terms of the 
variable m, since 2 now has another meaning. Thus we have m? — 44m + 
442 = 0, which has the double root m = 2; the solution to the differential 
equation is y = c,e** + c,xe***. Applying the boundary conditions and 
simplifying, we obtain 


c,=0 (1 +24) +¢,(22+ 24) =0 


It now follows that c, = 0 and either c, = 0 or 2 = —1. The choice 
c, = 0 results in the trivial solution y = 0; the choice 2 = —1 results in the 
nontrivial solution y = c,xe?*, c, arbitrary. Thus, the boundary-value 
problem has the eigenvalue 2 = —1 and the eigenfunction y = c,xe”*. 


iL x<2 


Solved Problem 15.3 Find a Fourier sine series for f(x) = ae on 
x 


(0,3). 


Using Equation 15.13 with L = 3, we obtain 


3 
C, = 3] Feosin a 


o7 NIX pe NIX 
= = | (0)sin — dx + = | (2)sin—dx 
3 5 3 3 5 3 
+ 3 al 4 2Qnx 
=O+ cos — cos cosnz 
3L na 3 j.0 na 3 
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Thus Equation 15.12 becomes 


4 2n1n nnx 
= 5 1)” {si 
F(x) 2 cos ( "Js 3 


2 
Furthermore, cos 3 = 5 cos a , COS ='1e.); 


Hence, 


4/1. am 3 .. 2m 2. 32x 
f(x)= sin sin +—sin as 
m2 3 4 3 3. 3 


Since f(x) is piecewise smooth on [0,3] and continuous everywhere in 
(0,3) except at x = 2, it follows from Theorem 15.3 that this equality is 
valid everywhere in (0,3) except at x = 2. 


Appendix 
LAPLACE 


TRANSFORMS 


Ax) F(s)=L4f0)} 

1. 1 — (s > 0) 
s 
1 

2, x — (s>0) 
Ss 

3 > GE 12d Sa (s > 0) 

s 


svn? (s > 0) 


5. vx 


Vas’? (s>0) 


OOG)---Qr-Wt nara 


6 G1) 7 
(s > 0) 
7 e (s >a) 
s-a 
8 “.— (s>0) 
sin ax ae 
s’+a 
Ss 

9, cos ax ae (s > 0) 

sta 
if 

10. sinh ax : (s > |al) 

s 
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ae 
Sx) F(s) = £{flx)} 
11. cosh ax 5 : 7) COG > lal) 
s°-a 
2as 
i s>0 
12 xX sin ax Gary ( ) 
so-a 
13. X COS ax s>0 
(s° + Ee ( ) 
n—1)! 
14, | oe” GH 1,2. 2.9 : : (s >a) 
s-a 
a 
1S; e”* sin ax Pe ee (s > b) 
s—b 
16. e”™ cos ax Gnae (s > b) 
2 
1. sin ax — ax cos ax | 5 = ay 9) 
(s° +a") 
18. | ose ! 
a 1+ as 
19 a “ 1) 
, i s(s —a) 
1 
20. 1- —x/a —S 
s(l+as) 
1 J te 1 
: a’ (1+as)? 
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Ax) F(s) = L4fo)} 
‘3 e® —e® ] 
. a=) (s—a)(s -b) 
53 exla etl 1 
. a—b (1+ as)(1+ 5s) 
S 
24. (1 + axje™ Gone 
i -x/a s 
25. —la—-x)je- “Say 7 OR 
a ( ) (i+ as)? 
ax —be® AY 
26. se ae = a SS 
a-b (s—a)(s —b) 
ae x! — be */4 Ss 
ab(a —b) (1+ as)(1+ bs) 
| peer: 1 
28. —(e™ -l-ax <a 
ae 3(5—a) 
i 2a? 
29. sin’ ax aa 
s(s° +4a*) 
a2 
30. sinh? ax —"- 
s(s° —4a*) 
cosh. eae 
1 43. ND a’ 
a V2 ax ax 44 a4 
sinh—= cos—= ae 
J2. V2 
‘ ax as 
32, sin —= sinh —= 
v2 v2 si +a‘ 


APPENDIX: Laplace Transforms 127 


x) F(s) = £4f@)} 
ax ax 
cos —= sinh —= + 
: v2 V2 as 
“a ay ax ax 4 4 
sin —— cosh — s'ta 
V2 Af 
ae Lees ee Ce (ee, 
3 
34. cos + cosh s 
ae: aa 
3 
35: J cinch ax — sin ax) fe 
2 nes 
2 
36. Gach ax — COS ax) as 
2 st_at 
2 
37. Sih ax + sin ax) as 
ms Cee 
1 @ 
38. — (cosh ax + cos ax) 
2 gage 
2 
i : 2a°s 
39. sin ax sinh ax ; 
s4+4a‘ 
2 2 
40. cos ax sinh ax a(s” —2a”) 
si 44a‘ 
2 2 
41. sin ax cosh ax ] +2a°) 
s*+4a‘* 
ie 
42. cos ax cosh ax 
s+ +44aqt 4 
1. 2 
43. |  —(sin ax + ax cos ax) ; a a 
2 (s° +a‘) 
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— area 
Ax) F(s)=L£6f)} 
3 
44. cos ax — sin ax 2 : 2 
2 (s +07) 
1 ; male Sa 
45. | —(axcoshax — sinh ax) 2 2980 
2 (s -a ) 
i: tie, as 
: a ax (s? -a2)y? 
2 
47. Gana oreddh ax) _ 2\2 
| 2 (s° -a‘*) 
3 
48 cosh ax + “sinh ax ay 
2 (s° -a‘) 
go,|  a8inbx —bsin ax e 
. Bee (s* +a°)(s? +b’) 
és cose SCosae — 
pia © (s? +a°)(s? +57) 
Al asin ax — b sin bx ca 
Ai ah? (s? +a7)(s? +57) 
- a* cos ax —b* cos bx s° 
; a? —b? (s? +a°)(s? +b?) 
? b sinh ax — a sinh bx ab 
53. eae (s? —a?)(s? —b?) 
; cosh ax — cosh bx 2. 
a gap (s° —a°)(s* —b*) 
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Ax) F(s) = L4fQ)} 
a asinh ax — b sinh bx s? 
* a2 —b* (s? —a*)(s? —b’) 
56 a’ cosh ax—b? cosh bx 3 
a? —b? (2 er SB) 
2 
57. ree a 
a s“(s° +a‘) 
est 
1 2 
58. —sinh ax — x ; . 5 
a s“(s* -a*) 
ax. a’ 
59. 1-cosax -—sin ax aa es 
2 s(s° +a“) 
ax. a’ 
60. 1—cosh ax + — sinh ax ares 
2 s(s* -a‘) 
| 24/2 
b? cosax —a* cosbx ab 
61. 1+ 5 ren 5 
a’ —b* S(s° +a°)(s° +b*) 
2 2 
a ee oe cosh bx azh2 
. a —b s(s? —a*\(s° =6*) 
63 1 (3-a*x’)sinax — ae 
8 3ax cos ax (s?4a7y 
aa ] a°s 
64. —|sin ax — ax cos ax —— 
8 (s° +a?) 
1 (l+a’x*)sin ax 7 a>s? 
65. — ae 
8 ax COS ax (s° +a‘) 
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es 


Ax) F(s) = L4fe)} 
66 1 (34+a’x*)sinh ax — a 
8 3ax cosh ax (ary 
3 
x a's 
67. — (ax cosh ax — sinh ax ——_— 
3! ) (s? ayy? 
a a ax cosh ax — as? 
; 8 (1—a*x”)sinh ax (s? — a’)? 
1 —x/ayn 1 
: —(1- 
i al ae S(as + 1)(as + 2)---(as +n) 
; ssinb+acosb 
70. sin (ax + 5) Tage a 
s“+a 
scosb—asinb 
71. cos (ax + 5) - © ei ae 
s“+a 
ax3 
cos - 5 
72 ax ax/2 2 3a 
_ axy3 sta 
43 sin 
2 
7 14+ 2ax Sta 
vinx svs 
1 
74, e/a 
sta 
AG 1 (e* ser) 
: xf s-a-vVs—b 
1 LS og /s 
16 cos2Vax —e" 
| vs 
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L I) F(s) = Lf} 
I 1 a/s 
v7, —=— cosh 2Vax Shug 
| ve aR 
1 
78. sin 2Jax —3/2 —a/s 
Van oes 
1 
79. sinh 2Vax -3/2 a/s 
Van Es Ne, 
i -a/s 
Jo(2Vax) 26 
s 
1 -als | 
Vx/aJ,(2Vax) se" 
AY 
sep (einer ay (2Vax) wee 
ae (p>9) | 
1 
83. Jn (x 
ue aes 
is 
$? Ay ag 
84. Ji(x) a ee 
s? +1 
(Vs? +1-5)? 
85. JAx) (p>) Ws tis)" 
Vs? +] 
1 
2a)’ T(p+— 
86. | x°J, (ax) [>> -4) (2a)°T(p 5 
: Vn(s? +q7)Pra/?) 
gPo 1 
87. >0 mailer 
T(p) (p > 0) a 
4" nl! n-(1/2) 1 
88. ————_— x 
(2n)!Vr tals 
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J) F(s) = £4fx)} 
go, | xe ga 0 i 
‘ e > 
T(p) eae (s +a)? 
ax dees 
90. l-e In S-a 
x S 
bx ax = 
91. e” -e In S-a 
x s-—b 
2. Sta 
92. — sinh ax In 
x S-a 
2 2 
93. Een in2 rile 
x 5° 
2 2 
94. S (cos bx — cos ax) In ical 
x 5° 4B? 
95. iat arctan — 
=| = : 
2. 2as 
96. —sin ax cos bx arctan 5 5 
x s°-a° +b 
a Pacer en? 
97. sin |ax| € ie Toeeae 
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Adams-Bashforth-Moulton 
method, 109, 111-112 

Amplitude, 53 

Analytic functions, 86—87 

Archimedes’ principle, 51 


Bernoulli equations, 6, 14 

Bessel functions, 99-100 

Boundary-value problems, 5, 
115-120 

Buoyancy problems, 51—52 


Characteristic equations, 34-35, 
73 
Circular frequency, 53 
Classifying solutions, 52—53 
Coefficients. See constant coeffi- 
cients; undetermined coeffi- 
cients; variable coefficients 
Completing the square, 58 
Constant coefficients, 30 
Laplace transforms solutions, 
66-67 
linear equations, 79-83 
Constant matrix, 70 
Convolutions, 59-60 
Critically damped motion, 53 


Damped motion, 53 

Decay problems, 20-21 
Defining properties, 10-11 
Denominators, 58—59 
Derivatives, 65-66, 72 


Differential equations. See also 
linear differential equations 
applications, 20-26 
classifying, 5-7 
definitions, 2—5 
first-order, 8—14, 20-26, 105-— 
112 
general solutions, 4—5 
notation, 3 
order, 3 
particular solutions, 4 
solutions, 3—5, 8-14 
Differential forms, 5—6 
Dilution problems, 23-25 
Direction fields, 105 
Dummy index, 101 


Eigenfunctions, 117, 119-120 
Eigenvalue problems, 117 
Eigenvalues, 117 
Electrical circuit problems, 49— 
51 
Electrical circuits, 25-26 
Equations. See also differential; 
linear differential 
Bernoulli, 6, 14 
characteristic, 34—35, 73 
exact, 7, 10-12 
homogeneous, 6—7, 9-10, 30, 
87-88 
indicial, 91 
nonhomogeneous, 32, 39-44, 
88-89 
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Equations (cont.) 
n-th order, 34—35 
ordinary differential, 2 
partial differential, 2—3 
second-order, 34, 35—36 
separable, 7 
Equilibrium position, 47 
Euler’s method, 106, 108, 111 
Exact equations, 7, 10-12 


Falling body problems, 21—22 
First-order differential equations, 
8-14 
applications, 20—26 
numerical methods, 105-112 
First-order systems, 79-83, 110 
Forced motion, 53 
Fourier series, 120-121 
Free motion, 52—53 
Froebenius, method of, 90-91 
Functions 
Bessel, 99-100 
gamma, 98-99 


Gamma functions, 98—99 
General solutions, 4—5, 32 
Growth problems, 20-21 


Homogeneous equations, 6—7, 9- 
10, 30, 87-88. See also linear 
differential equations 

Hooke’s law, 48—49 


Identity matrix, 72 

Independent variables, 57 

Indicial equation, 91 

Initial-value problems, 5, 9, 44, 
82, 89 

Integral, 72 


Integrating factors, 11-12 

Inverse Laplace transforms, 57— 
58 

Isoclines, 105 


Kirchhoff’s law, 50 


Laplace transforms, 56-62, 124- 
132 
inverse, 57-58 
solutions by, 65-67 
Laws 
Hooke’s, 48—49 
Kirchhoff’s, 50-51 
Newton’s of cooling, 21 
Newton’s second of motion, 22 
Linear differential equations, 6, 
12-14 
first-order, 105-112 
homogeneous, 30, 34-37, 39-— 
44 
nonhomogeneous, 39-44 
reduction to first-order system, 
79-83 
second-order, 47—53, 86, 115— 
117 
solutions, 34—37, 79-83 
theory of solutions, 29—32 
Linearly dependent solutions, 31 
Linearly independent solutions, 
31 
Linear systems, 67 
Line elements, 105 


Matrices, 70—75 

Matrix addition, 71 

Matrix exponential, 73-75 
Matrix multiplication, 71 
Matrix solution methods, 79—83 


Method of Frobenius, 90-91 
Method of undetermined coeffi- 
cients, 40—42 
Methods 
Adams-Bashforth-Moulton, 
109, 111-112 
Euler’s, 106, 108, 111 
Frobenius, 90-91 
modified Euler’s, 108 
numerical for first-order equa- 
tions, 105—112 
predictor-corrector, 108 
Runge-Kutta, 108-109, 111 
undetermined coefficients, 40— 
42 
Modified Euler’s method, 108 
Motion, 52—53 


Natural frequency, 53 

Natural length, 49 

Newton’s law of cooling, 21 
Newton’s second law of motion, 
22 

Nonhomogeneous equations, 32, 
39-44, 88-89 

Notation, 3 

nth-order equations, 34—36 
Numerators, 59 

Numerical instability, 107 
Numerical methods, 105-112 


Ordinary differential equations, 2 
Ordinary points, 87 

Orthogonal trajectories, 26 
Oscillatory damped motion, 53 
Overdamped motion, 53 


Partial differential equation, 2—3 
Partial fractions, 58—59 
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Particular solutions, 4 
Period, 53 
Phase angle, 53 
Power series, 86—92 
Power series method, 87—88 
Predictor-corrector method, 108 
Problems 
boundary-value, 5, 115-120 
buoyancy, 51-52 
decay, 20-21 
dilution, 23-25 
eigenvalue, 117 
electrical circuit, 49-51 
falling body, 21-22 
growth, 20-21 
initial-value, 44, 82 
spring, 47-49 
Sturm-Liouville, 117-118 
temperature, 21 


Recurrence formula, 88 

Regular singular points, 89-90 

Runge-Kutta method, 108-109, 
111 


Scalar multiplication, 71 
Second-order equations, 34, 35- 
36 
Second-order linear differential 
equations, 47-53, 86, 115- 
117 
Separable equations, 7 
Simple harmonic motion, 53 
Solutions, 3—5 
classifying, 52—53 
first-order differential equa- 
tions, 8-14 
initial-value problems, 82 
Laplace transforms, 65-67 
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Solutions (cont.) 
linear differential, 34—37, 79— 
83 
linearly dependent, 31 
linearly independent, 31 
linear systems, 67 
matrix methods, 79-83 
no initial conditions, 82—83 
particular, 4 
power series, 86—92 
theory of, 29-32 
Spring problems, 47—49 
Square matrix, 70 
Standard forms, 5—6 
Steady-state motion, 53 
Step sizes, 107 
Sturm-Liouville problems, 117-— 
118 


Temperature problems, 21 
Transient motion, 53 
Translations, 61-62 


Undamped motion, 53 
Underdamped motion, 53 
Undetermined coefficients, 40—42 
Unit step function, 60-61 


Variable coefficients, 30, 86 
Variation of parameters, 42—44 


Vectors, 70 


Wronskian, the, 31—32 


